CONSTRUCTION OF THE PAULI-VILLARS-REGULATED 
DIRAC VACUUM IN ELECTROMAGNETIC FIELDS 
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Abstract. Using the Pauli-Villars regularization and arguments from 
convex analysis, we construct solutions to the classical time-independent 
Maxwell equations in Dirac's vacuum, in the presence of small external 
electromagnetic sources. The vacuum is not an empty space, but rather 
a quantum fluctuating medium which behaves as a nonlinear polariz- 
able material. Its behavior is described by a Dirac equation involving 
infinitely many particles. The quantum corrections to the usual Maxwell 
equations are nonlinear and nonlocal. Even if photons are described by a 
purely classical electromagnetic field, the resulting vacuum polarization 
coincides to first order with that of full Quantum Electrodynamics. 
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1. Introduction 

In classical Physics, a time-independent external density of charge p ex t 
and a charge current j C xt induce a static electromagnetic field (E ext = 
— VFcxt, -Bcxt = curl^ext); which solves Maxwell's equations in Coulomb 
gauge: 

-Ay cxt = 47repext, 
< -AA ext = 47rej'ext, (1-1) 
div A ext = 0, 
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where e is the elementary charge. It is convenient to gather the electro- 
static and magnetic potentials in a unique vector A ext = (V ex t, A cxt ) called 
the four-potential, which we will do in the whole paper. The electromag- 
netic potential A cxt solving (1.1) is the unique critical point of the (time- 
independent) Maxwell Lagrangian action functional 

j^Pcxtjcxt^) = — [ (\VV(x)\ 2 - | cmlA(x)\ 2 ) dx 
8tt J R 3 V / 

-el p ext (x)V(x)dx + e j cxt (x) ■ A(x) dx, (1.2) 
Jr 3 Jr 3 

which is strictly convex with respect to V and strictly concave with respect 
to A. In particular we can obtain A ext by a min-max procedure: 

3f p extJext (A ext ) = mmmaxJz? PextJext (V, A) = maxmmJz? PextJext (V, A) 
V A A V 

where the constraint div A = is always assumed. 

The situation is much more complicated in Dirac's vacuum. It has been 
known for a long time that, in reality, the vacuum is not an empty space, but 
rather a quantum fluctuating medium which behaves as a nonlinear polariz- 
able material [14, 13, 27, 36, 21]. In this medium, virtual electron-positron 
pairs induce a polarization in response to external fields. The resulting elec- 
tromagnetic field which is observed in experiments has to take into account 
the vacuum polarization effects. The corresponding four-potential A* solves 
coupled nonlinear Maxwell equations of the form 

' - = 4tt e (/9 V ac(eA*) + p ext ) , 

< -AA* = 4vre (j vac (eA,) + j cxt ), (1.3) 
div A* = 0. 

Here p va c(eA*) and jvac(eA*) are respectively the charge density and the 
charge current induced in the vacuum. As we shall explain, they are non- 
linear and nonlocal functions of eA*. 

The Dirac vacuum is described by Quantum Field Theory, that is, by a 
second-quantized fermion field. The charge and current densities / o vac (eA*) 
and j V ac(eA*) are obtained by minimizing the energy of this field in the 
presence of the given potential eA*. In this model the interaction between 
the Dirac particles is mediated by the classical electromagnetic field which 
accounts for photons. This approach to vacuum polarization is usual in the 
Physics literature (see, e.g., [36, 20]). 

The main idea behind the nonlinear Maxwell equations (1.3) is that the 
vacuum behaves as a nonlinear medium, and the form of the equation is 
reminiscent of nonlinear optics. The nonlinear effects are in practice rather 
small since e has a small physical value, but they become important in the 
presence of strong external sources, which can produce electron-positron 
pairs in the vacuum. Already in 1933, Dirac computed in [9] the first or- 
der term obtained by expanding /> vac (eA») in powers of e. The nonlinear 
equations (1.3) was then studied by Euler, Heisenberg, Kockel and Weis- 
skopf among others [14, 13, 27, 46]. In a celebrated paper, Schwinger [36] 
used (1.3) (and a time-dependent version of it) to derive the probability of 
pair creation by tunneling in a strong electrostatic field. For more recent 
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works on the subject, see the references in [20]. Several ongoing experiments 
aim at detecting some nonlinear effects of the vacuum in the laboratory [6] . 

Like for the usual Maxwell equations (1.1), the nonlinear equations (1.3) 
in Dirac's vacuum arise from an effective Lagrangian action, which now 
includes the vacuum energy 



Here J-" vac (eA) is the ground state energy of Dirac's vacuum in the potential 
eA. The densities of the vacuum are then defined by 

ep vac (A) := — T vac (eA) and ej VSLC {A) := - — T vac (eA). (1.5) 

We note that the vacuum correction — J r vac (eA) to Maxwell's Lagrangian has 
been computed to first order in the semi-classical approximation in [28, 36]. 

It is not so easy to provide a rigorous definition of the vacuum energy 
J r vac (eA). It is well-known that this quantity is divergent in the high energy 
regime and an ultraviolet regularization has to be imposed. In this paper we 
use the famous Pauli-Villars regularization method [32] to properly define 
the vacuum energy J 7 vac (eA) (Theorem 2.1 below). Then we are able to 
state our main result (Theorem 2.4 below), which gives the existence of 
a critical point of the effective Lagrangian action (1.4), when the external 
sources p cxt and j cxt are not too large. As a corollary, we obtain solutions 
to the nonlinear Maxwell equations (1.3). 

This article is the continuation of several works dealing with the Hartree- 
Fock approximation of Quantum ElectroDynamics (QED), some of them 
in collaboration with Solovej, [22, 23, 25, 24, 18, 19], and which originated 
from a seminal paper of Chaix and Iracane [4]. There, only the purely 
electrostatic case was considered. To our knowledge, the present work is the 
first dealing with electromagnetic fields in interaction with Dirac's vacuum. 

Acknowledgements. M.L. and E.S. acknowledge support from the French 
Ministry of Research (Grant ANR-10-0101). M.L. acknowledges support 
from the European Research Council under the European Community's 
Seventh Framework Programme (FP7/2007-2013 Grant Agreement MNIQS 
258023). We are also grateful to the referees for their careful reading of the 
manuscript. 



2.1. Elementary properties of electromagnetic Dirac operators. Be- 
fore entering the main subject of this article, we recall some elementary 
spectral properties of the Dirac operator in the presence of electromagnetic 
fields [43, Chap. 4]. 

We work in a system of units such that the speed of light and Planck's 
reduced constant are both set to one, c = H = 1. We introduce the Dirac op- 
erator with mass m, elementary charge e and electromagnetic four-potential 
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A = (V,A), 

D m ,eA := a ■ ( - «V - eA(x)) + eV(x) + m/3 (2.1) 

which is an operator acting on L (R 3 ,C ). Here the four Dirac matrices 
a. = (ai, CK2, 0:3) and /3 are equal to 

- - u ?) - " - (» a x) • 

the Pauli matrices 0*1, er 2 and 0-3 being defined by 

<n:=(; J), <r 2 :=(° and <r 3 := (j ^ . 

The spectrum of the free Dirac operator is not semi-bounded [43], 
c(-Dm,o) = ( — 00, — m] U [m, 00). 

As we will recall below, the unbounded negative spectrum of D m< o led Dirac 
to postulate the existence of the positron, and to assume that the vacuum is a 
much more complicated object than expected. The mathematical difficulties 
arising from the negative spectrum are reviewed for instance in [12]. 

In our setting, the natural space for the four-potential A = (V, A) is the 
Coulomb-gauge homogeneous Sobolev space 

Hi v (R 3 ) := [A = (V,A) e L 6 (R 3 ,M 4 ) : 

divA = and F = (-VV.curlA) G L 2 (M 3 ,M 6 )}, (2.2) 
endowed with its norm 

ll^ll^i iv (R3) := II w IIl2(]r3) + II curl ^lli 2 (R 3 ) = II- f IIl2(r3)- ( 2 - 3 ) 

Here and everywhere, the equation div^4 = is understood in the sense 
of distributions. The requirement F G L 2 (IR 3 ) simply says that the elec- 
trostatic field E = —W and the magnetic field B = curl^4 have a finite 
energy, 

/ \E 2 \ + \B\ 2 < 00. 

Lemma 2.1 (Elementary spectral properties of D^a)- Let m > 0. 

(i) Any four-potential A 6 i/j; (K 3 ) is D m $-compact. The operator D m A 
is self-adjoint on H 1 ^ 3 ) and its essential spectrum is 

a css (D„ hA ) = (-00, -m] U [m,oo). 

(ii) The eigenvalues of D m< A in (—in, m) are Lipschitz functions of A in 
the norm ^ R3 y 

(Hi) There exists a universal constant C such that, if 

ll^lljfi.^MS) < vVm, (2.4) 
for some number rj < 1/C, then 

<T(D m>A ) n (-m(l - Cn), (1 - Cn)m) = 0. 
(iv) Finally, ifV = 0, then o-(D m> A) H (— m,m) = 0. 
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Proof. Recall the Kato-Seiler-Simon inequality (see [37] and [39, Thm. 4.1]) 

\/p>2, ||/(*M-iV)|| e^-^H/IUMU (2-5) 

(2tt)p 

where & p is the usual Schatten class [39]. Applying (2.5) with p = 6 together 
with the Sobolev inequality, we obtain 



V- 1 



D 



m,0 



and, similarly, 



ct • A- 



D 



C C 

< -— |M| < _||w|L 2 , 



C 1 1 ,n c „ , 

r~ — ~~/ = \\ \\l 6 — ~ 7 = CU1 " L 2 ' 

6 6 v m v m 



m,0 

where we have used that div^4 = 0. Since all the operators in ©6 are 
compact, statements (i) and (n) follow from usual perturbation theory [29, 
35]. Concerning (in), we notice that 

D m , A (D^o)- 1 = (l+(V-a-A) 

so that, under condition (2.4), 

\D mA \ > (l-C^)|D m , |. 
Statement (in) then follows from (i), whereas (iv) is [43, Thm 7.1]. □ 

2.2. The Pauli-Villars-regulated vacuum energy. 

2.2.1. Derivation. We are now ready to define the energy of the vacuum, 
using Quantum Field Theory. We consider a fermionic second-quantized 
field, placed in a given electromagnetic potential A. Later, the potential A 
which describes light and external sources will be optimized. In this section 
it is kept fixed and we shall look for the ground state energy of the Dirac 
field in the given A. Our second-quantized field only interact through the 
potential A, there is no instantaneous interaction between the fermions. 
The Hamiltonian of the field reads 



H 



eA 



■=\ I fa*(x)D m!eA y(x)-y(x)D m , eA y*(x))dx, (2.6) 
2 Jr3 V / 

where *(x) is the second-quantized field operator which annihilates an elec- 
tron at x and satisfies the anti-commutation relation 

+ = 2S a , u S(x - y). (2.7) 

Here 1 < a, v < 4 are the spin variables and ^(x) a is an operator- valued 
distribution. The Hamiltonian M. eA formally acts on the fermionic Fock 
space T = C © Ar>1 /\^ L 2 (M 3 , C 4 ). The proper interpretation of the 
expression in parenthesis in (2.6) is 

**(x)D m:eA V(x) - y(x)D m:eA y*(x) 

4 

:= ^(x),(D m , eA )^(x) u -^(x) tl (D m , eA )^*(x) u . 

/U,I/=1 
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This choice is made to impose charge-conjugation invariance [25], 

OT^ -1 = HT eA 

where ^ is the charge-conjugation operator in Fock space [43]. 

For any fixed A, the Hamiltonian (2.6) is a quadratic polynomial in the 
creation and annihilation operators ^f*(x) and *&(x). It is therefore very 
well understood. The expectation value in any state in Fock space can be 
expressed as 

(U eA ) = tr D m<eA - (2.8) 

where 7 is the one-particle density matrix of the chosen state, namely 

l(x,y)a,u = (**(z) CT *(y)„). 

The subtraction of half the identity comes from charge-conjugation invari- 
ance. The details of this calculation can be found in [25]. Because electrons 
are fermions, it is known that 7 must satisfy the constraint < 7 < 1 on 
L 2 (M 3 , C 4 ), which is called the Pauli principle. Conversely, any operator 7 
such that < 7 < 1, arises (at least formally) from one state in Fock space. 
Since the energy only depends on the quantum state of the electrons through 
the operator 7, we can refrain from formulating our model in Fock space and 
only use the simpler operator 7 and the corresponding energy (2.8). 

We note that the energy (2.8) is gauge invariant. Namely, it does not 
change if we replace A by A + Vx and 7 by e tex je~ tex for any function %. 

Remark 2.1. A preferred state among those having 7 as one-particle den- 
sity matrix is the unique associated quasi- free state [2] , also called (general- 
ized) Hartree-Fock state. So the main simplification with the model of this 
paper is that, when the electromagnetic field is purely classical, the ground 
state of the Hamiltonian is always a quasi-free state. Hartree-Fock theory 
is exact here, it is not an approximation. This simplification does not occur 
when the photon field is quantized. 

We are interested in finding the ground state of the vacuum, which cor- 
responds to minimizing (2.8) with respect to 7. For atoms and molecules, 
we would impose a charge constraint of the form 

tr( 7 -i)=iV. 

In this paper we restrict ourselves to the vacuum case for simplicity, and we 
thus do not have any other constraint than < 7 < 1. The formal minimizer 
of the energy (2.8) is the negative spectral projector 

7 = l(-oo,0)(An,eA)- 

The interpretation is that the polarized vacuum consists of particles filling 
all the negative energies of the Dirac operator D my eA-, in accordance with 
the original ideas of Dirac [7, 8, 9]. 1 

^For atoms and molecules, the vacuum projector l(-oo,o)(' ' ' ) nas to be replaced by a 
spectral projector of the form l(_oo,fx)( - • ' )i f° r some chemical potential /i which is chosen 
to ensure the correct number N of electrons in the gap (more precisely the correct total 
charge of the system). Except from this change of chemical potential, the equations take 
exactly the same form. 
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The corresponding ground state energy is 

<!^<i tr Dm > eA ( 7 ~ 2~) = ~ \ tr I Dm ' eA I • ( 2 - 9 ) 

Of course this energy is infinite, except if our model is settled in a box with 
an ultraviolet cut-off [25]. 

Here and everywhere in the paper, the absolute value of an operator is 
defined by the functional calculus 

|C| := VC*C. 

It is in general not a scalar operator, that is, it may still depend on the spin. 
In the special case of D m fi, it does not depend on the spin, however, since 
it is the scalar pseudo-differential operator 

|An,o| = V-A + m 2 . 

In order to give a clear mathematical meaning to (2.9), we argue as follows. 
First, we can subtract the (infinite) energy of the free Dirac sea and define 
the relative energy as 

F rc \{eA) := -tr(|l> m ,o| " \D m ,eA\). (2.10) 

Since we have removed an (infinite) constant, we formally do not change the 
variational problem in which we are interested, hence we also do not change 
the resulting equations. 

Unfortunately, the functional (2.10) is not yet well-defined, because the 
model is known to have ultraviolet divergences. Indeed, the operator |-D m ,o|— 
\D„ h eA\ is not trace-class when A / 0, which is reminiscent of the fact that 
the difference of the two corresponding negative projectors is never Hilbert- 
Schmidt [31]. This can be seen by formally expanding the trace in a power 
series of eA. The first order term vanishes and the second order term is 
computed in Section 4.2 below. It is infinite if no high energy cut-off is 
imposed. 

In order to remove these divergences, an ultraviolet cut-off has to be 
imposed. The choice of this regularization is extremely important. When 
the trace of |-D m ,o| — |-Dm,eA| is expanded as a power series of eA, several 
terms which look divergent actually vanish because of gauge invariance. In 
addition to obvious physical motivations, it is necessary to keep the gauge 
symmetry for this reason. Some simple choices in the spirit of what we have 
done in the purely electrostatic case (see, e.g. [18] for two different choices) 
would not work here, because of their lack of gauge symmetry. 

In 1949, Pauli and Villars [32] have proposed a very clever way to regular- 
ize QED, while keeping the appropriate invariances. It is this technique that 
we will use in this paper (but there are other choices, like the famous dimen- 
sional regularization [30] ) . The Pauli- Villars method consists in introducing 
J fictitious particles in the model, with very high masses mi, mj playing 
the role of ultraviolet cut-offs. 2 These particles have no physical significance 
and their role is only to regularize the model at high energies. Because of 



Since in our units c = h = 1, m has the dimension of a momentum. 
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their large masses, they do not participate much in the low energy regime 
where everyday Physics takes place. 

In our language, the Pauli-Villars method consists in introducing the fol- 
lowing energy functional 



(2.11) 



see Remark 2.2 below for some details. Here tuq = m and cq = 1, and the 
coefficients Cj and rrij are chosen such that 




J J 

E^ = E c ; m i = - ( 2 - 12 ) 

3=0 3=0 

It is well-known in the Physics literature [32, 20, 3] that only two auxiliary 
fields are necessary to fulfill these conditions, hence we shall take J = 2 in 
the rest of the paper. In this case, the condition (2.12) is equivalent to 



ci 



ml 



777*2 — 171*1 ,,tj 2 — 

We will always assume that tuq < ra\ < 777-2, which implies that c\ < and 
c 2 > 0. 

The role of the constraint (2.12) is to remove the worst (linear) ultraviolet 
divergences. In the limit 777,1,7772 — > 00, the regularization does not prevent 
a logarithmic divergence, which is best understood in terms of the averaged 
ultraviolet cut-off A defined as 



and c 2 = 



2 2 

m\ — TTIq 



777o 



777 



2 ' 



(2.13) 




(2.14) 



The value of A does not determine mi and 7772 uniquely. In practice, the 
latter are chosen as functions of A such that c\ and C2 remain bounded when 
A goes to infinity. 

That the model is still logarithmically divergent in the averaged cut-off 
A can again be seen by looking at the second order term in the expansion, 
given by (2.20) and (2.22) below. Removing this last divergence requires a 
renormalization of the elementary charge e. This can be done following the 
method that we used in the purely electrostatic case with a sharp cut-off 
in [19], but it is not our goal in this paper. 

Remark 2.2. In our language the fictitious particles of the Pauli-Villars 
scheme are described by density matrices 7^, with 70 = 7 and the divergent 
energy of the vacuum is chosen in the form 

J 



tTD m . ieA (ij - 1/2) 



3=0 

instead of (2.8). When optimizing the energy subject to the Pauli principles 
< 7j < 1, one has to minimize over the matrices jj such that Cj > 
and maximize over those such that Cj < 0. Adding the infinite constant 
Y. J j=o c j tr|-Dmj,o|/2 gives Formula (2.11). 
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2.2.2. Rigorous definition. Our main result below says that the energy func- 
tional A i y J r pv(A) is well defined for a general four-potential A = (V, A) 
in the energy space i^ iv (M 3 ) (and which therefore satisfies the Coulomb 
gauge condition div A = 0) . 

Theorem 2.1 (Proper definition of Tpy). Assume that Cj and mj satisfy 

2 2 

Co = 1, ?7l2 > 777-1 > 777-0 > and Cj = Cj777| = 0. (2-15) 

j=0 j=0 

(i) Let 

1 2 

T A ■= - 5>(|A^,o| - \D mj , A \). (2.16) 

3=0 

For any A € L^R^R 4 ) n ij^R 3 ), the operator tr C 4 T A is trace-class on 
L 2 (R 3 ,C). In particular, J 7 py(A) is well-defined in this case, by 



Jpv(A) = tr(tr C 4TA). (2.17) 



(n,) TTte functional Tpy can be uniquely extended to a continuous mapping 
onH^R 3 ). 

(in) Let A G i^ iv (R 3 ). We have 



Fpv(A)=T 2 (F) + K(A), (2.18) 



where 



where F := (E,B), with E = — W and B = curl A The functional 1Z is 
continuous on H\ iv (R 3 ) and satisfies 

|KWI „((gM)|| F || l2 + (gN) ||F|| ^, (2 . 19) 

for a universal constant K . 

(iv) The functional Ti is the non-negative and bounded quadratic form on 
l A ) given by 

MF) = ^J R3 M(k)(\B(k)\ 2 -\E(k)\ 2 )dk, (2.20) 

2 2 r 1 

M(k) := J^Cj u(l - it) log [m] + u(l - u)\k\ 2 ) du. (2.21) 

The function M is positive and satisfies the uniform estimate 

< M(k) < M(0) = 2log ( A \ (2.22) 

3tt 

where A was defined previously in (2.14). 

Let us emphasize the presence of the C 4 -trace in statement (i) about 
the trace-class property of tr^T^. We do not believe that the operator 
is trace-class without taking first the C 4 -trace, except when V = 0. If we 
are allowed to take more fictitious particles by increasing the numbers of 
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auxiliary masses, it is possible to obtain a trace-class operator under the 
additional conditions 

Yl ''■> '" ' = Yl C 'J " r j = "• 

3 j 

At least four auxiliary masses are then necessary. The terms which are not 
trace-class when only two fictitious particles are used, actually do not con- 
tribute to the final value of the energy functional J-'py (their trace formally 
vanishes). For this reason, we have found more convenient to first take the 
C 4 -trace (which is enough to discard the problematic terms) and limit our- 
selves to two fictitious particles, as is usually done in the Physics literature. 
This suffices to provide a clear meaning to the energy. 

The function M describes the linear response of the Dirac sea. It is well- 
known in the Physics literature [20, Eq. (5.39)]. We will see below that 

lim f^A _ M{k) \ = m := \jf f l * dz . (2.23) 

A^ooV 3vr v ') v 1 4vr J q 1 + |fc|2(l - z 2 )/4 V ; 

The function in the right-hand side of (2.23) was first computed by Ser- 
ber [38] and Uehling [44]. The same function U already appeared in our 
previous works dealing with pure electrostatic potentials [26, 23, 19]. This 
is a consequence of the gauge and relativistic invariances of full QED. 

Our proof of Theorem 2.1 in Sections 3 and 4 below, consists in expanding 
the energy J r py(A) in powers of the four-potential A. We use the resolvent 
expansion at a high but fixed order and therefore our main result is valid for 
all A's, not only for small ones. All the odd order terms vanish (by charge- 
conjugation invariance). The second order term is given by the explicit 
formula (2.20) and it is responsible of the logarithmic ultraviolet divergences. 
It will be important for our existence proof that this term be strictly convex 
in A and strictly concave in V. We also have to deal with the fourth order 
term in some detail. The latter was computed in the Physics literature 
in [28] and our task will be to estimate it. The higher order terms are 
then bounded in a rather crude way, following techniques of [22]. The main 
difficulty in our work is to verify that the Pauli-Villars conditions (2.12) 
induce the appropriate cancellations in the few first order terms, and to 
estimate them using the L 2 -norm of the electromagnetic fields and nothing 
else. 

In spite of its widespread use in quantum electrodynamics, the Pauli- 
Villars scheme [32] has not attracted a lot of attention on the mathematical 
side so far (see [16, 40, 41, 42, 15] for some previous results). Theorem 2.1 
seems to be among the first in this direction. 

2.2.3. Differentiability. After having properly defined the functional J-py, 
we need some of its differentiability properties, in order to be able to define 
the charge and current densities of the vacuum by (1.5). In this direction, 
we can prove the following 

Theorem 2.2 (Differentiability of Tpy). Assume that Cj and rrij satisfy 
conditions (2.15). 
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(i) Let A G i7^ iv (M 3 ) be such that is not an eigenvalue of the operators 
D mjt A for j = 0, 1, 2. Then the functional JFpy is C°° in a neighborhood of 
A. 

(ii) The first derivative of Fpy is given by 

<dF PV (A),(t>,a)> = / <(pA,-M,),(t>,a)) K4) (2.24) 

for all (t>, o) G i^ iv (R 3 ), where the density p A and the current j a are defined 
as 

Pa{x) := [tr C 4Q A ](x,x) and j A (x) := [tr C 4 a Q A ] (x, x), (2.25) 
and with Qa refering to the kernel of the operator 

2 

Qa ■= 1 (-00,0) (An.,, a)- 

3=0 

The operators tr C 4 Qa and tr C 4 ci^Qa for k = 1,2,3 are locally trace-class 
on I? (M 3 , C 4 ) , and pa and ja are well-defined functions in L\ oc (M?) n C, 
where C is the Coulomb space 

\m\ 2 



:= {/ : M 3 -»• C : / 



dk < oo 



(2.26) 



/rs |Af 

(m) There exists a universal constant rj > swc/t i/tai i/te second derivative 
of Fpy satisfies the estimate 

< 2K ( y 1 -^- 



dlFp V (A) - i- 



-M 
M 



3=0 3 



IAII 2 

div 



(2.27) 



for all A such that \\A\\^i ^ ffi3 ^ < r/^mo = n^fm. 



div 



Our estimate (2.27) means more precisely that 

<A',d 2 Jp V (A) A'} - ±- J^M(k)(\B'(k)\ 2 - \E'(k)\ 2 ) dk 



<2k( y^i 



i=o 



I . i|2 

1 11 "div 



lA'l 



^div 



when A is small enough in H^ iy ( 

As a consequence of Lemma 2.1 and Theorem 2.2, we obtain 

Corollary 2.1 (Regularity in a neighborhood of 0). There exists a positive 
radius rj such that the functional J-pv is C°° on the ball B(rf) := {A G 

#divC* 3 ) 



|A||^i ( R 3) < ry^/mo}. On this ball, the differential d-Fpv is 
given by (2.24), whereas d 2 J 7 pv satisfies estimate (2.27). 

Proof. We fix rj such that 

Cr,< 1, 

where C is the constant in statement (ra) of Lemma 2.1. For this choice, 
given any four-potential A in the ball B(rj), is not an eigenvalue of each 
of the operators D mjjA . Corollary 2.1 then follows from Theorem 2.2. □ 
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2.3. Solutions to the Maxwell equations in Dirac's vacuum. In this 
section, we explain how to use Theorems 2.1 and 2.2 in order to get the 
stability of the free Dirac vacuum, and to construct solutions to the nonlinear 
Maxwell equations. The proofs of Theorems 2.1 and 2.2 will be detailed 
later. 

Let e > be the (bare) charge of the electron. Assume that cq = 1, 
and that Cj and rrij satisfy (2.12). We work under the condition that e < e 
for some fixed constant e (e is not allowed to be too large, but it can be 
arbitrarily small). All our constants will depend on e, but not on e. Note 
that e is dimensionless here because we have already set the speed of light 
equal to 1. 

Using Theorem 2.1, we can properly define the effective electromagnetic 
Lagrangian action by 



for all A 6 i^ iv (IR 3 ), the Coulomb-gauge homogeneous Sobolev space. Our 
purpose will be to construct a critical point of this Lagrangian, which will 
in the end solve the nonlinear equations (1.3). 

Remark 2.3. In this paper we have considered a second-quantized Dirac 
field which only interacts with a classical electromagnetic field. There is an- 
other way to arrive at exactly the same Lagrangian action (2.28), which is 
closer to our previous works [22, 23, 25, 24, 18, 19]. We start with Coulomb- 
gauge QED with quantized transverse photons. Then, we restrict our at- 
tention to special states in Fock space of the form £1 = S7hf <8> ^Coh, where 
Ohf is a n electronic (generalized) Hartree-Fock state characterized by its 
one-particle density matrix < 7 < 1, and Slcjoh is a coherent state for 
the photons, characterized by its magnetic potential A(x) (a given classical 
magnetic potential in R 3 ). We thereby get a Hartree-Fock model coupled to 
a classical magnetic field. Because of the instantaneous part of the Coulomb 
interaction, the model contains an exchange term. When this term is ne- 
glected, one obtains the exact same theory as in this paper. In relativistic 
density functional theory [11], the exchange term is approximated by a local 
function of p 7 _i/2 and j 7 -i/2 only. 

2.3.1. Stability of the free Dirac vacuum. In the vaccum case p ext = Jcxt = 0, 
we have the obvious solution A = 0. The following theorem says that is 
stable in the sense that it is a saddle point of the effective Lagrangian action, 
with the same Morse index as for the classical Maxwell Lagrangian action. 
This can be interpreted as the stability of the free vacuum under its own 
electromagnetic excitations. 

Theorem 2.3 (Stability of the free Dirac vacuum). Assume that Cj and rrij 
satisfy (2.15). The four-potential A = is a saddle point of J^py It is the 
unique solution to the min-max problem 



j^° v (0,0) = min _~2p V (^0)= max _^ Y (0,A), (2.29) 





curl A\\ L 2 < 



e 
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or, equivalently, 

jzf° v (0,0)= min sup Sf$y(V,A) 

livv|| i2 <r^s || curlA || L2< i^i 

max inf _JSf# v (V,A), 

||curlA|| L2 <I^0 || W || i2 <I^a 



(2.30) 



for some positive radius r which only depends on Y^j=o \ c j\( m o/ m j) an d e 
(the largest possible value of e). 

The result is a direct consequence of the properties of the functional Ti 
defined in (2.20), as well as of the regularity properties of J-pv- 

As we have seen we can let the cut-off A go to oo (which implies that 
mo/nij —7- for j = 1,2), while keeping c\ and C2 bounded. We therefore 
see that the radius r of the ball of stability of the free vacuum does not go to 
zero in the limit A — > oo if the bare parameters e and mo are kept bounded. 

For A = 0, the electrostatic stability of the free Dirac vacuum was pointed 
out first by Chaix, Iracane and Lions [4, 5] and proved later in full generality 
in [1, 22, 23]. It is possible to include the exchange term and even establish 
the global stability of the free Dirac vacuum [22, 23, 24]. Dealing with 
magnetic fields is more complicated and, so far, we are only able to prove 
local stability, using the Pauli-Villars regularization. Because of lack of 
gauge symmetry, it is not clear whether the free Dirac sea is still stable 
under magnetic excitations when a sharp ultraviolet cut-off is used. 

Proof. We choose r > such that 

r<ri/V2 and 2k( V ^] m (r 2 + 2m r 4 ) < — ^, (2.31) 

m jJ 8vre 2 

where K is the constant appearing in (2.19), and where r\ is the constant 
in statement (Hi) of Theorem 2.2. We recall that e < e. Consider now any 
A such that || W||x,2 < r^/mo/e and || curl A\\ L 2 < r^/mo/e (which implies 
H^llifi. (r 3 ) = II-^IIl 2 ^ \/2m r/e). By (2.19), we have 

< 2k( ^^jmo(r 2 + 2m r 4 )e 2 ||F||2 2 



- 8jr \\ \\L* 



(2.32) 



Using Formula (2.20) for J^, we get 

^p°v(U, 0) > |^ f M(k)\E(k)\ 2 dk > 0, 
with equality if and only V = 0, since M > 0. Similarly, 

=2 



^pv(0,A) <~ f M(k)\B(k)\ 2 dk < 0, 
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with equality if and only A = 0. Thus we have shown (2.29). The equiva- 
lence between (2.29) and (2.30) is a classical fact of convex analysis, see [10, 
Prop. 1.2, Chap. VI]. 

Finally, since we can deduce from (2.27) that 



V z — ' rrij I 

V 7=0 ■> 7 



2 

L 2 



'3 = 

for e||A||^i , ffi3 N < ry/2m,Q < rj^/rriQ, we deduce that Jz?p V is strictly convex 
with respect to V and strictly concave with respect to A, provided that r 
satisfies the additional condition 

AKmJy^\ 2 <^. (2.33) 
This implies uniqueness of the saddle point by [10, Prop 1.5, Chap. VI]. □ 



2.3.2. Solution with external sources. We now come back to our initial prob- 
lem and consider an external density p ext and an external current j cxt . It 
will be convenient to express our result in terms of the size of 

V'cxt := ep C xt * rr and A cxt := ej ext * t^-, 
\x\ \x\ 

which are the corresponding potentials when the vacuum does not react. 
We look for a saddle point of the Lagrangian action jf^™ t, -' oxt defined above 
in (2.28). 

Theorem 2.4 (Nonlinear Maxwell equations in small external sources). 
Assume that Cj and mj satisfy (2.15). Let r be the same constant as in 
Theorem 2.3. 

(i) For any 

e||AKt||fli v(R 3) < -*g— , (2-34) 

there exists a unique solution A* = (V*,A*) G iJ^ iv (IR 3 ) to the min-max 
problem 



(A*) = min _ j^f>^ (y, A,) 
l|w|| L2 <^o 

max A), 



(2.35) 



r s /rriQ 



||curlA|| i2 < ^ 

or, equivalently, to 

^Poxt,ioxt( A ^ = min sup jgf^^(A) 

Iiw|| i2 <^g5 ||curlA||i2< i^o 

max inf _ J^ xt ,icxt (A). 



(2.36) 
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(ii) The four-potential A* is a solution to the nonlinear equations 

\ -AV* = Aire (p eAt + p ext ) , (9 „ 7 a 

\ -A^=47re(j e ^+j cxt ), 

where p e A, and j e A„ refer to the charge and current densities of the Pauli- 
Villars-regulated vacuum 

2 

Q* = ^2 Cj l(-oo,0) {D mj ,eA,), (2.38) 
3=0 

defined in Theorem 2.2. 

Equations (2.37) and (2.38) are well known in the Physics literature (see, 
e.g., [11, Eq. (62)-(64)]). Solutions have been rigorously constructed in the 
previous works [22, 23, 25], with a sharp ultraviolet cut-off, but in the purely 
electrostatic case A ext = A* = 0. In this special case it is possible to obtain 
the polarized vacuum as a global minimizer. The method of [22, 23, 25] does 
not seem to be applicable with magnetic fields, however. To our knowledge, 
Theorem 2.4 is the first result dealing with optimized magnetic fields in 
interaction with Dirac's vacuum. 

The proof of Theorem 2.4 is based on tools of convex analysis, using that 
^Pextjext l oca l saddle point geometry by Theorem 2.1. 

Proof. Let us define the balls 

By(r) := {V G L 6 (M 3 ,M) : e||W|| L 2 < r^}, 

and 

B A (r) := {A G L 6 (M 3 ,M 3 ) : e|| curlA|| L 2 < r^rri^}. 

As we have already shown in the proof of Theorem 2.3, when r satisfies 
condition (2.31), the function A i— >■ «ifp V (A) is strictly convex with respect 
to V and strictly concave with respect to A on By(r) x B\(r). We deduce 
that 

A ' y J^v(A) + e / j cxt • A - Pcxt V 

satisfies the same property. 

We now assume that the external field A cxt G By(er) x B^(er) for some 
e < 1/3 to be chosen later, and we look for a saddle point in £>v(3er) x 
£>A(3£r). Since J^y'' 3 ' 6 * 4 is strongly continuous on £> v (3er) x £>A(3er) by 
Theorem 2.1, a classical result from convex analysis implies that jgP^J^ 1 
possesses at least one saddle point A* = (V^A*) G £>v(3er) x B\(Ser), 
solving 

j^ext,jext/ A j = min £>p™^(V,AA= max j£ftf J " rt (V;, A). 
v uee v (3 £ r) FV A€,B A (3£r) ^ v 

See for instance [10, Prop. 2.1, Chap. VI]. Uniqueness follows from the 
strict concavity and convexity, by [10, Prop 1.5, Chap. VI]. 
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It only remains to verify that A* does not lie on the boundary of By(3er) x 
Bx(3er). Similarly as in (2.32), we first compute 

|^Pv(eA') " ^(eF')| < k(± M) (e 4 ||F'||l 2 + ^ll^ll^) 

S=o m ^ V m ° / (2.39) 



(3e) 



2 



< ^ — —\\F' 
~ 8vr 



/||2 



2 



for all A' £ By(3er) x £ A (3er), when r satisfies (2.31) and e < 1/3. We 
obtain with £> C xt := curlvlext 

<- c 2jr 2(0jB )_J_ / 5 2 + -3-/ B^-B+^f 2E 2 + B 

87T J R 3 47T J R 3 87T J R 3 

When ^4 belongs to the boundary of BA(3er), we obtain 

||£ ext || L2 <l||5|| L2 and \\E\\ L2 <\\B\\ L2 

and therefore 

J^Pext Jext (y A) _ ^Pext Jext (y Q) < — e 2 J r 2 (0, B) - 1 ~ 2 / 3 ~ 3 ( 3g ) 2 /" #2 

Choosing e = 1/9, the right-hand side is < — e 2 J r 2(0, B) < since B / 
0. So we have shown that when A belongs to the boundary of By(3er), 
£>P^>j^(y,A) < J^ xt ' icxt (F,0). Since A* maximizes A t-> Jz?^ xt ' jcxt (K, A), 
it cannot have an energy smaller than ^fpy t '"' ext (K lt , 0) and we conclude that 

e 2 / | curl.4*! 2 < (3e) 2 r 2 m . 

Similarly, we can show that K does not belong to the boundary of By{3sr). 

The unique saddle point A* = (V*,A*) being in the interior of the set 
£>v(3er) x £>A(3er), the derivative of j^P^' 1 ^ must vanish at this point. 
Using the value (2.24) of the derivative of Tpy computed in Theorem 2.2, 
we find that 

f - AT4 = 47T e (p eA , + Pcxt) 

I -A A* = 4vr e (j eAt + jcxt) 

where p ej 4 t and j e A, are given by (2.25) in Theorem 2.2. This is exactly the 
self-consistent equation (2.37). □ 

The rest of the paper is devoted to the proofs of Theorems 2.1 and 2.2. 
Our strategy is as follows. First, in Section 3, we show that the functional 
J-pv is well-defined for four-potentials A with an appropriate decay in x- 
space (the integrability of A on IR 3 is enough). Then, we compute things 
more precisely in Section 4, and we exhibit the cancellations which show 
that this functional can be uniquely extended by continuity to H^ iy ( 
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3. The Pauli-Villars functional for integrable potentials 
The purpose of this section is to prove that the operator 

1 2 

tr C 4 T A := -^Cjtr C 4 (\D mjfi \ - \D mjjA \), (3.1) 

3=0 

is trace-class, when the four-potential A := (V, A) decays sufficiently fast. 
The proof relies on an expansion of J-py(A) with respect to the four- 
potential A using the resolvent formula, but for which we actually do not 
need that A is small. Our precise statement is the following 

Proposition 3.1 (tr C 4 Ta is in Assume that Cj and rtij satisfy con- 

ditions (2.15). Then, the operator tr C 4 Ta is trace-class whenever A € 

l x (]R 3 ,ir 4 ) nff 1 ^ 3 ^ 4 ). 

Remark 3.1. For this result, it is not important that div^4 = 0, hence we 
do not require that A G H\ iv ( 



The rest of this section is devoted to the proof of Proposition 3.1. 
Proof. Our starting point is the integral formula 

^2 



X 



1 f x , 1 f ( iuj iuj \ , . 
= - / — ^doj = — / 2 + )doj. (3.2) 

JR X Z + UJ Z 27T J]R V X + IUJ X — IUJ J 



T 2 +oj 2 



When T is a self-adjoint operator on L 2 (IR 3 ,R 4 ), with domain D(T), it 
follows from (3.2) using standard functional calculus (see e.g. [34]), that the 
absolute value |T| of T is given by 

T = — / ( 2 — — + - )dw. (3.3 

1 1 2vr J R V T + ioj T - iuj) y ' 

Let us remark that this integral is convergent when seen as an operator from 
D(T 2 ) to the ambient Hilbert space. In particular, 

< min jl,u; 2 ||^ 2 || I )( r 2)^ i 2( R 3 C 4) [• 

Since the domains of D^. q and A are both equal to H 2 (M. 3 , C 4 ), we 
deduce that we can write 

1 r 2 

m 1 / \ - / IUJ IUJ 

T A = — / > Cj( 

4 W K ^ \D mj}A + iuj D mjA -iuj 

iuj iuj \ , 
- — + — duj 

^rrij ,0 

+ iuj D mi p-iu)J 

on i/ 2 (M 3 ,C 4 ). Here and everywhere else it is not a problem if D mjt A has 
as an eigenvalue. The operator D mjt A + iu is invertible for to / 0, and 
(iuj) (D mj a + iu)^ 1 stays uniformly bounded in the limit uj — > 0. 

In order to establish Proposition 3.1, we will prove that the C 4 -trace of 
the integral in the right-hand side of (3.4) defines a trace-class operator 
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according to the inequality 



IOJ IUJ 



L g cjtrc4 (- 

lit) it,) \ 

duj < oo, (3.5) 



. () -II,. A + I In A ~ iu 



IUJ IUJ 

D mj:0 + ioj D mjfi -iuj, 



©i 



which we can establish when A G L 1 (R 3 ,1R 4 ) n ^(R 3 , M 4 ). This will com- 
plete the proof of Proposition 3.1. 

As a consequence, our task reduces to derive estimates in Schatten spaces 
on the integrand operator 



K{u, A) := Cj tr c 



IOJ IUJ 



D mi , A + D m A - ioJ 



IOJ IOJ 

+ 



Dm^Q + iuj D m0 -ioj 



)- 



which we can integrate with respect to oj. To this end, we use the resolvent 
expansion 

ioj ioj %oj / , \ 1 
— — = — (a ■ A — V) — , 

D mj ,A + t^J D mj j0 + ioj D mj > A + ioj D mj j0 + ioj 

iterated six times and obtain 

5 . , 

IUJ IUJ v- ioj ((a A V) Y 

D mj ,A + iu D mj>0 + iuj ^D m3i0 +iwv D mjfi + iuJ 

+ 7T^ ({cc-A-V) — ^ ) 6 , 

D mj ,A + iOJ V D mj:0 + iojJ 

(3.6) 

together with the similar expression for the term with — ioj instead of +ioj. 
Again, we insist on the fact that this expansion makes perfect sense for 
oj 7^ 0, even if the spectrum of D m a contains 0. This allows us to write 



K(oj, A) = Y^ tr C 4 (^R n (oj, A) + R n {-oj, A)) 

+ tr C 4 (R' 6 (oj,A) + R' 6 (-oj,A)), 



n=l 



(3.7) 



with 

2 



and 



Rn(u,A) :=r&- — ((a-A-V)- — Y , (3.8) 



R' 6 (uj,A):=y Cj - -f( a .A-V)- -) . (3.9) 
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Our purpose is to prove that 

5 



[ ( J2 || tr c4 (Rn(w, A) + R n (-u, A)) || 

■/K\ n= l 

+ || tr C 4 (R^(cu, A) + R'si-cu, A)) \\ & ^ du> < oo. (3.10) 

Estimate on the sixth order term. We first estimate the sixth order 
term R' e (to, A) in (3.7) which is the simplest one. The C 4 -trace is not going 
to be helpful for us here. First we use the inequality 



<1, 



which, in particular, takes care of the possibility of having in the spectrum 
ofD mj A- Combining with Holder's inequality in Schatten spaces, ||AB|| 6i < 

le„ W B lev' we obtain 



IbJ 



D m] ,A + ibJ 



Si 



< 



(a ■ A - V) 



1 



D m -o + iu 



(3.11) 



e 6 



We next use the Kato-Seiler-Simon inequality (2.5), similarly as in the proof 
of Lemma 2.1, which gives us 

1 



Vp>3, 
where 



(a-A-V) 



6, 



Ir, :-- 



1 



< (JpHmf + w 2 )2p 5 ||A|| LP , 



r 2 dr 



2vr 2 7o (l + r 2)l' 
For p = 6, we can use the Sobolev inequalities 

\\V\\ L e < S\\VV\\ L 2 and ||A|| £8 < S\\VA\\ L 2, 
to obtain an estimate in terms of the gradient VA by 



(3.12) 



( a -A-V)—- 



+ iui 



Inserting in (3.11), we have 



® 6 (m 2 +w 2 )4 



, ii va Hl- 



so that 



/ 



j=o (m 2 + w 2 )2 
\R' 6 (u;,A)\\ 6l ^ < S G lJ^ l ^)\\VA\\% f - 



. j= 0"-J/ ^(1+W 2 )2 

The term with +ioj replaced by — ioj is treated similarly. 



(3.13) 



(3.14) 
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Estimate on the fifth order term. The method that we have used for 
the sixth order term of (3.6) can be applied in a similar fashion to the fifth 
order term, leading to the estimate 

2 

. I r ■ I \ _ / 1 1 rlf , ? 

(3.15) 

None of these estimates use simplifications coming from the C 4 ~trace. The 
latter is only useful for lower order terms. 



2 

/ \\R 5 (±u,,A)\\ 6l du<lJ^^)\\A\\ 5 L5 / 77 



Estimate on the fourth order term. For the other terms in (3.6), we 
need more precise estimates based on conditions (2.13) satisfied by the co- 
efficients Cj and the masses rrij. We start by considering the fourth order 
term, for which we use the identity Co + c\ + C2 = to write 

2 4 



j=0 k=0 m °' u 



( 



ICO 



D 



)(( a .A-V) - 1 . ) 



4-fc 



Next we use that 
ioj 



IOJ 



mo — mj 



(3- 



IUJ 



D mj fl + ioj An ,o + ioj 
'uj\ 



< (rrij — mo) 



TTIq + 0J A 



since rrij > mo, and we argue as before, using this time A £ L l 
obtain 

2 



\\R 4 (±c,A)\\ ei 



< 



5I4\oj\ 



(m§ + w 2 )2 j= o 



{rrij-mo) ||A||| 4 , 



(3.16) 



. We 



(3.17) 



hence 



/ ||i? 4 (± W ,A)|| ei ^<5/ 4 ||A||l 4 ^| Cj |^_^ / (3.18) 

Notice again that we have not used the C 4 -trace in our estimate of the 
fourth order term. 



Estimate on the first order term. In order to deal with the first, second 
and third order terms, we need to use more cancellations. We start by 
considering the first order term for which we can write 



IUJ 



2lo 2 



(a ■ A - V) 
(a ■ A - V) 



IUJ 



D 



m,-,0 



+ iu D r . 



+ 



(a ■ A - V) 
(a ■ A - V) 



Inserting 



D 



D mA n ± iu 



rrij ,0 T 



1 



D 



rrij ,0 — 



2cj 2 



(3.19) 
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we obtain 



( a .A-y) - — — - - - - (cL-A-vy 



D mjfi + iuj y D m . fi + iuj D mjfi -iuj y >D mjfi -iuj 

2u 2 , .1 

where the notation {T\,T2} refers to the anti-commutator operator 

{T U T 2 } :=T X T 2 + T 2 T X . 

At this stage, we recall that 

(a -X) (a-Y) = X Y + i(X x Y) • S, (3.20) 

for all (X,Y) £ (K 3 ) 2 . In this formula, X x Y is the cross product of 
the vectors X and Y, whereas the notation S = (Si, S2, S3) refers to the 
matrices 

>:: ",)■ < 321 > 

As a consequence, we obtain 

{a-p,a-A} = {p, A} R3 + i(p x A + A x p) • S = {p, A} R3 + B ■ S, 

where {•, -} R 3 is a notation for 

{S,T} R3 :=S-T + T-S 

and where p = — iV (a simplifying notation that will be used henceforth). 
Since f3atk + = 0, we deduce that 

{at - A — V, D m ., } = { P ,A- V<x} R3 + B ■ S - 2m j F/3. (3.22) 

This finally gives us 

22i(w,j4) + fli(-w,j4) = 2w 2 (^ l5l +^1,2), (3.23) 

where 

2 1 1 

^1,1 : =E c i^ T—^({P^ A - Va }^+ B -' E )— I T—^> ( 3 - 24 ) 

z —' J p 2 + mj + u! 2 \ l JK J p 2 + mj + uj 2 

]=0 J 1 3 



and 



Ki, 2 :=-2Vc 3 m JT - ^_^/3 (3.25) 

^ p 2 + mj + u 2 p 2 + mj + lo 2 
3=0 3 3 

Concerning the operator 1Z\^, the last step consists in using identities (2.12) 
and the two expansions 

1 1 m\-m 2 , 

+ 



p 2 +m 2 -+uj 2 p 2 + m\ + oj 2 (p 2 + m\ + oj 2 )(p 2 + m 2 - + uj 2 ) 



1 mi - mj , 

+ 7 2^ 2 ' 2,2 ( 3 ' 26 ) 



p 2 + m,Q + cj 2 (p 2 +TOq + w 2 ) 



(mo - m 2 ) 2 



(p 2 + m 2 + oj 2 )(p 2 + m\ + oj 2 ) 2 



22 P. GRAVEJAT, C. HAINZL, M. LEWIN, AND E. SERE 

This gives 

2 



n 1A = J2 c i( m o - m i) 2 ( 

o— n V 



3=0 

; ({p,A-Va} R3 +B--£) 



(p 2 + tuq + uj 2 ) 2 ' ! ~ (p 2 + m\ + cj 2 )(p 2 + m 2 + uj 2 ) 

+ (pi+m^+u^+ml+uj 2 ) 2 ^ A ~ Va}R3 + B ' T ' ) p 2 + m 2 + uj 2 

+ 2 I \ I 2 (fa ^ " + B ■ S) 



P 2 + 77T.Q + W 2 ' "' (p 2 + TO 2 + 0J 2 )(p 2 + Wq + w 2 ) 2 

We now use the fact that A 6 L 1 (IR 3 , M 4 ), B = i(pxi + ixp), as well as 
the Kato-Seiler-Simon inequality (2.5) to get 

2 j 
||fti,i|| 6l < 18g |c 3 -|(mg - m 2 ) 2 K + 7 ^ 2)2 || A|| L1 . (3.27) 

The analysis of the operator lZ± t 2 is more involved. Under conditions (2.12), 
we are not able to prove that 1Z\ t 2 is trace-class. However we can compute 
first the C 4 -trace before taking the operator trace. We obtain 

tr C 4 fti l2 = 0, (3.28) 

since tr C 4 (3 = 0. 

Remark 3.2. By this argument, we do not prove that IZi^ is trace-class. 
Under the additional conditions 



the operator IZip becomes a trace-class operator, and its trace is equal to 
0. This strategy however requires to introduce additional fictitious particles 
in our model. Introducing more fictitious particles in order to justify the 
computation of a term which is anyway does not seem very reasonable 
from a physical point of view. This explains why we prefer here to first take 
the C 4 -trace. 

As a consequence, we can conclude our estimate of the first order term 
by combining (3.27) and (3.28) in order to obtain 



f ti C 4(R 1 (cu,A)+R 1 (-u,A)) 

JR 



6i 



< 367 7 V \cA— — A L i / 5^. 3.29 

j m o h 1+w 2 2 
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The second and third order terms are treated following the same method, 
except that the algebra is a little more tedious. We start by writing that 

2 



R n (oj, A) + R n (-cu, A) = 2oj 2 £ Cj £ ( 1 , (a ■ A - V)) 

j=0 k=o KIJm iV w J 

D 2 m . fi + u 2 V >D mjfl + iu) 

We next expand as before using (3.19). 



k 

x 



Estimate on the second order term. For the second order term, we are 
left with 



R 2 (cu, A)+R 2 (-uj, A) = -2oj 2 ^ c 

3=0 



■3 X 



^ 2 // . T T \ 1 \2 A/ D m .Q 



x I —2 ^1 (a 



A - v )7^—,) -£ 



X (a - A - V))*^^( (a • „ - V)^^^" 



which may also be written as 

2 



R 2 (u, A) + i? 2 (- W) A) = -2u 2 £ cj ( 2 1 (a ■ A - y) 2 1 

J=0 V % ,0 + w ^mj ,0 + c 



(3.30) 
Inserting 

(a. - A - V) 2 = \ A\ 2 + V 2 - 2 a • A V, 
and (3.22) into (3.30), we are led to 

R 2 (u,A) + R 2 (-gj,A) = -2uj 2 (K 2tl + K 2<2 ), (3.31) 

where 

2 



1 / 1 
^2,1 :=y^Cj-5 s ? ( (\A\ 2 + V 2 -2a- Av)-„ 2 , 

^— ' p + mf + w z \ ^ ' P + mi + v 

j — 3 j 

- Am2 3 ( V 2 , \ , 2 ) 2 " ( ({?> A " FQ }r3 + B ■ S) ^ 2 ) ^ , 



3 

and 



x ^ 1 

^2,2 :=2^c j m j ^ 2 + ^ 2 

<f^/3 2 ^ V 2 ,({p,A-Fa} R3 + fl-s) 1 2 |. 

[ p l + m| + ui z \ k J M / + + J 



p z + + uj 2 X 
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The proof that 72.2,1 is trace-class is similar to the first order case, us- 
ing (3.26). The final estimate is 



||rc2,i|| 6l <£M s/ 7 



3=0 

X 



m" — nin 



\B\\ T 2+h-+ 



m - — rrir 



(4m|||y||^ 2 + ||B||* a + 8(m 2 - m 2 )\\ A\\\ 2 ) ) . 



Since 

tr C 4 7e 2 ,2 = 0, 
as for the first-order term, our final estimate is 



(3.32) 



/ ||tr C 4 (R 2 {oj,A)+R 2 (-uj,A)) 

JR 11 



6i 



3=0 



1 1 

ml. — rrin 



m 



u r dco 

A hl\\ B \\Ll L (1 +W 2)2 



+ h 



m; — m. 



o 



x (4m|||^ 2 + |K 2 + 8K 2 -^)||A||y / -^rY (3-33) 

Estimate on the third order term. Similar computations give for the 
third order term 

2 2u 2 



R 3 (u, A) + R 3 (-u, A) = J2 <H ( p 2 + Z+A ^ A " Va ^ 3 + B '^ 

3=o i 

~ *™,V0) . * 2 ) 3 - - 2 -4t A (\A\ 2 + V 2 -2 a -AV) 
' p 2 + m 2 + oj 2 J p 2 + m 2 + uj 2 [ v ' 



-> 2 + m 2 + uj 2 

p 2 + mf + w 2 k J M y p 2 + mf + uj 2 



(3.34) 



Using once again (3.26), we deduce 



f tr C 4 (Rs(uj, A) + i? 3 (-w, A)) duj<KS2\ Cj \ ( 
Jr ©i Ji 



da; 
(l+w 2 )f 



x ( m^sllVHia + /sllVH^IIBII^ + ^^||5||| 2 ||F|| i4 + ^ll||B||| 2 



m ■ 



+ m j 7 7 ||A|| i 3||F||| 3 + { ^tfi\\B\\ 2 L2 \\A\\ Li + A|| A||| 3 ||F|| L a 



ra- 



mi 



+ (m ■ - m ) 



^r4llA||| 4 ||S|| L2 + ^ 



for some universal constant K. 



Combining with (3.14), (3.15), (3.18), (3.29) and (3.33), we obtain (3.10), 
provided that A is in L^R 3 ,^ 3 ) n ^(R^R 3 ). This concludes the proof of 
Proposition 3.1. □ 
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4. Estimates involving the field energy 
In Proposition 3.1 above we have shown that the operator 

1 2 

tr c4 T A := - tr C 4 (\D mjfi \ - \D mj>A \), 

3=0 

is trace-class when A decays fast enough. More precisely, in the proof of 
Proposition 3.1, we have written 

T A = V T n (A) + I* (A) :=-?- V / (i? n (w, A) + i? n (-^, A)) d W 

47r n^i^ (4.1) 

+ -/-/ (R' 6 {cu,A) + R' 6 {-u;,A))cbj, 
4vr J R 

with i? n and R' 6 given by (3.8) and (3.9), and we have proved that the op- 
erators tr C 4 T n (A) and tr^4Tg(A) are trace-class. However our estimates 
involve non gauge-invariant quantities (some IP norms of A) and they re- 
quire that A decays fast enough at infinity. 

In this section, we establish better bounds on these different terms. We are 
interested in having estimates which only involve the field F = (— VV, curl A) 
through the norms ||W||£2 and || curl A||^2. Our simple estimate (3.14) on 
the sixth order only depends on the field F. But we will also need to know 
that the sixth order is continuous, which will require some more work. For 
the other terms, we have to get the exact cancellations. 

With these estimates at hand, it will be easy to show that J^y can be 
uniquely extended to a continuous function on the Coulomb-gauge homoge- 
neous Sobolev space ijj iv (M 3 ), as stated in Theorem 2.1, and which we do 
in the next section. 

Remark 4.1. In the estimates of the previous section, it was not important 
that div ^4 = 0. We have to use this property now. 

4.1. The odd orders vanish. The following lemma says that the trace of 
the odd order operators tr C 4 Ti(A), tr C 4 T%(A) and tr C 4 T§{A) vanish. This 
well-known consequence of the charge-conjugation invariance is sometimes 
called Furry's theorem, see [17] and [20, Sec. 4.1]. 

Lemma 4.1 (The odd orders vanish). For A G i^ iv (M 3 ) n L 1 (R 3 ,M 4 ) and 
n = 1, 3, 5, we have 

tr ( tr C 4 T n {Aj) = i- J tr ( tr C 4 (R n (cj, A) + R n (-u, A))) duj = 0. (4.2) 

Proof. Let Cip := ifta,*^ be the (anti-unitary) charge-conjugation operator. 
Since CD mj , C _1 = -D mjJQ , we have 

C(D m . fi ±iu)~ 1 C- 1 = -(A^oiiw)" 1 . 
Similarly, since A and V are real-valued, we can write 

Ccx- AC' 1 = a- A and CVC~ 1 = V, 

so that 

C Ru(±oj, A) C" 1 = {-l) n R n {±Lu, A). (4.3) 
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At this stage, we can compute 

tr c4 (CTC- 1 ) =tr c4 T, (4.4) 
for any operator T on L 2 (IR 3 ,C 4 ). Here, T refers to the operator defined as 

T{f) :=TQ). 

When tr C 4 T is trace-class, so is the operator tr C 4 T, and its trace is equal 

to 

tr ( tr C 4 T) = tr(tr C 4 T). (4.5) 

As a consequence, the operator tr C 4(CTC~ 1 ) is trace-class, as soon as T is 
trace-class, and its trace is the complex conjugate of the trace of T. 

Finally, recall that we have established in the proof of Proposition 3.1 that 
the operators ti C i(R n {uj, A) + R n (—uj, A)) are trace-class for n = 1,3,5. 
Combining (4.3) with (4.4) and (4.5), we obtain 

tr(tr C 4 (R n (u,A) + R n (-oj,A))^ 

= (-l) n tr(tr C 4 (R n (u,A) + R n (-co,A))y 

We deduce that the quantity tr ( tr C 4(i? n (w, A) + R n {— w,A))) is purely 
imaginary when n is odd, so that the trace of tr C 4 T n (A) is purely imaginary. 
Since the operator tr C 4 T n (A) is self-adjoint, its trace is necessarily equal to 
0. This gives Formula (4.2). □ 

4.2. The second order term. We now compute exactly the second order 
term T 2 {A) appearing in the decomposition of Ta, assuming that A belongs 
to i7 1 (lR 3 ,M 4 ) and div A = 0. We will verify that it only depends on the 
electromagnetic fields E := — W and B := curl A. 

Lemma 4.2 (Formula for the second order term). For A £ i?^ iv (IR 3 ) n 
L 2 (R 3 ,M 4 ), we have 

tr(tr C 4T 2 (A)) = -L / M(k)(\B(k)\ 2 - \E(k)\ 2 ) dk := T 2 (F), (4.6) 

where M is the function defined in (2.21) and F = (E, B). 

The formula (2.21) for the function M which is nothing but the dielectric 
response of Dirac's vacuum is well-known in the physical literature. The 
following proof is inspired of the calculations in [20, p. 280-282]. 

Proof. In the course of the proof of Proposition 3.1, we have shown that 
the operator tr C 4T2(A) is trace-class when A G ^(IR 3 ,!^ 3 ) (see inequal- 
ity (3.33)). As a consequence, its trace is well-defined and given by 

tr(tr c4 T 2 (A)) = I (tr^2KA))(p,p)dp. (4.7) 

Here, tr^T^A) refers to the Fourier transform of the trace-class opera- 
tor tr C 4X2(A). Our convention for the Fourier transform of a trace-class 
operator T is the following 

T(p,q):=-^ [ T(x,y)e-^e^ydxdy. 
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In view of (3.31), the operator tr C 4 T 2 {A) is given by 
tr c .T 2 (A)(p,p) = - [ f> 2 ,\, 2 (-{\A\ 2 + V 2 )- 



+ ({P, -4} R 3^— 2~ 2~) + 5^- Bfe_ 2 _ ^ 2~ 2 Bk ~T~, 2 

V <- JI pHm,!^/ f-^ p z + m z - + uj z p z + 



+ {pfc,^}^- — — 2"{p fc ' y }^ — ~2~ — 2 

^— ' J p z + + ' p z + m z - + uj z 

k=l J ' J 

1 2 \ 

+ Am2 j( V — l 2—^) Wd^. 

J V p 2 + m z + w 2 / / 



(4.8) 



Since clival = 0, we deduce after a lengthy calculation that 

3 

tr(tr c4 T 2 (A)) =^T 2 , fe , (4.9) 

where 



d/c w 2 iiw dp 



Co - 



(4.11) 



Q J (p 2 + m 2 + w 2 ) 2 ((p — k) 2 + m 2 + w 2 ) 

x ((p • %)) (p • A(-k)) + (p 2 + m?)|y(fc)| 2 ) 
:=T 2 , 2 (.4)+T2,2(n 

and 

_ 1 f A fc 2 |l(A;)| 2 + (fe 2 -4p-A:)|y(A : )| 2 „ 2j J 

T 2 ,3:=— A Yl °i 12! \ \ vk<( TxTZ 2) 2^ dk u duj d P- 

47r 4 J R 7 ^ (p z + m z + uj 2 ) 2 ((p - k) 2 + m z + uj 2 ) 

(4.12) 

We next use the following Ward identities [45] 



PmPn dp 



J (p 2 + m 2 + o; 2 ) 2 ((p - k) 2 + m) + u 2 ) 

(4.13) 



X 3 S 



, (fcm ~ 1m){K ~ q n )dq 

' j ((« - *0 2 + ™? + w 2 ) 2 (g 2 + m 2 + w 2 ) ' 



for all (m,n) £ {1,2,3} 2 and all k £ M 3 . This equation is nothing else 
than a change of variables p = k — q, which makes perfect sense thanks 
to conditions (2.12) which guarantee the convergence of the integral. Its 
importance is well-known in the Physics literature, see, e.g., [33, Sec. 7.4]. 
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Since divA = 0, we infer that 

3 3 

X L § ° jPmdpn ((p 2 + m 2 + u J 2 )((p-k) 2 + m 2 + u J 2 )) ^ 
Integrating by parts, we are led to 

/2 ' 21 j 47r 4 7 M7 ^ Cj (p2 +m 2 +a; 2 )((p _ A . ) 2 +m 2 + a; 2 ) - ^ J 

Similarly, we can compute 



\V{k)\ 2 dkuj 2 dudp 



Q ' J (p 2 + m 2 + w 2 )((p - /c) 2 + m 2 + w 2 ) 

2 |F(A;)| 2 dfea; 4 dwdp 



vr 4 7 M 7 ^ (p 2 + m 2 + w 2 ) 2 ((p - /c) 2 + m 2 + w 2 ) ' 



Integrating by parts with respect to u, one can check that 

oj a dm dp 



E 



^ (p 2 + m 2 + w 2)2((p - t) 2 + " 2 + " 2 ) 



3 /" v—-v 



u; 2 dp 



Q ' 3 (p 2 + m 2 + 0J 2 )((p - k) 2 + m 2 + oj 2 ) ' 



so that 

1 /" A |V"(A;)| 2 dA;fa; 2 da;dp 

T 2 ,2(V) " 4^4 _/ M7 C i (p2 + m 2 + W 2)((p _ jfe)2 + m 2 + W 2) ' ( " 5 ^ 

On the other hand, since A and V are real-valued, we have 

|Ap(0) + F 2 (0) = - —-3- / (|A(fc)| 2 + |t>(fc)| 2 ) dfc, 
(27r)2 Jr3 

hence 

_ 1 r A |l(fc)| 2 + |y(fc)| 2 2 , „ , 

T 2 ,i = -— / V c, ' \ 1 ' \ ; ' 2 a; 2 dw dk dp. 
47r A 7 (p 2 +m 2 + u 2 ) 2 



r 1 f (2p-fc-|fc| 2 )(|l(fc)| 2 + |y(fc)| 2 ) 

72,1+72,2 = -^ / ZCjT^- 2~ — 2\2?7 2~^ — 2\ 03 dudkdp. 

4?r JWT? n (p 2 + m 2 + uj 2 ) 2 ((p - k) 2 + m 2 + oj 2 ) 



Combining with (4.14) and (4.15), we arrive at 

{2p.k_-\k\ 2 )(\A(k)\ 2 + \V(k)\ 2 ) 

In view of (4.9) and (4.12), this provides 

tr(tr c4 T 2 (A)) = / G(k)(\A(k)\ 2 - \V(k)\ 2 ) dk, (4.16) 
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where 



We next use the identity 



p • ku 2 doj dp 

Co 



Q 3 (p 2 + m 2 + uj 2 ) 2 ((p — k) 2 + m 2 + (j 2 ) 



-L = jf S 2 e - S(u a + (1- U)6 ) d< ^ u 

see [20, Chap. 5], to rewrite 
1 f 2 

G(k) = — j I y Cj uj 2 du> p ■ k dpx 

27r Jr 4 



i ' (4.17) 




J 



Using conditions (2.12), we can invoke Fubini's theorem to recombine the 
integrals in (4.17) as 

i t>l poo 2 

G{k) = ±.J J Y^ 3 z~ S{ > +{l ~ U)) s 2 dsudux 



I p-k e s(j> 2 -2(l-u)p-k) f j e sJ>Jl duj) dp. 

Jm. 3 \Jr J 



Since 



/>«" Vdu, = ^|, 

JR 2S2 

and 

/ p- ke- s{p2 - 2{l ~ u)p - k) dp=k- V( / e P-xsp 2 -Mi-u) P -k dp 

= (^) l (i-„)fcV( 1 -«) 2fc: 



x=0 

2 1,2 



we deduce that 



£2 /•! roo 2 



G[k) = ^j / Vc j e- s K 2 +«( 1 -«) fe2 ) s - 1 d S u(l-n)^. (4.18) 



i=o 

Integrating by parts, we now compute 

2 



poo 
J ,_ n 



c . e -K+«(i-«)^) 5 -i ds 

j=0 



/•CO _ 

/ V Cj log(s)e- s K 2 +«( 1 -«) fc2 ) (m 2 + «(1 - u)k 2 ) ds, 
Jo J=0 
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which is justified again thanks to conditions (2.12). Letting a = s(m 2 + 
u(l — u)k 2 ), we infer again from (2.12) that 

poo 2 2 

/ E^^ 1 ""' k) s- 1 ds= / y>i og (— — n — ui)^^ 

Jo f^ Jo f^ \m 2 + u(l -u)k 2 J 

2 

= ~y^ y Cj \og(m 2 + u(l - u)k 2 ). 

3=0 

Inserting into (4.18), we get 

k 2 f 1 k 2 
G ( k ) = ~-£~2 J 2Z c i u ( l - u ) log(w- +u(l - u)k 2 )du = —M(k). 

Combining with (4.16), we obtain Formula (4.6). □ 

We complete our analysis of the second order term by giving the main 
properties of the function M. 

Lemma 4.3 (Main properties of M). Assume that Cj and rrij satisfy (2.15). 
The function M given by (2.21) is well-defined and positive on M 5 , and 
satisfies 

< M(k) < M(0) = 2 -^-, 
where A is defined by (2.14). Moreover, 

2 1og(A) \k\ 2 f 1 z 2 -z 4 /3 



f 1 z l -z 

J TTJkW 



3vr w 4vr J I + \k\ 2 (l - z 2 )/A 

when mi — > oo and mi — > oo. 

Most of the above properties of M are well-known in the Physics litera- 
ture, see for instance [20, Sec. 5.2]. The positivity of M(k) for all k, which 
is crucial for our study of the nonlinear Lagrangian action, does not seem 
to have been remarked before. 

Proof. In view of (2.21), the function M is well-defined on ]R 3 . Concerning 
its positivity, we set 

2 

~Y. c i teg(m 2 + t), 

3=0 

for all t > 0. Using (2.13), we compute 

&( t ) = V Cj = ( m i ~ m o)( m i ~ m o) > o 
^m 2 +t (ml + t){m\ + t)(m\ + t) 



Since 3>(0) = -2 log A < and 



by (2.12), we deduce that 

-2 log A < $(t) < 0, 
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for all t > 0. Inserting into (2.21), we obtain (2.22). 
As for (2.23), we first write 

M(k) = - ^ jf * J> u(l - u) ( log(m 2 ) + log (l + " (1 ~ 2 " )fc2 )) d« 
2 log A 2 f 1 A . / w(l-u)/e 2 \ 

J — J 

When mi — > oo and r»2 — > oo, we infer that 
2 log A 



(fat. 



M(k) 



3vr 

Integrating by parts, we compute 



2 f 1 

-/ u(l - u) log (1 + u(l - u)k 2 ) du. 
77 Jo 



[ u(l-u) log (l + u(l - u)k 2 ) du = - [ 
Jo Jo 



'- 1 (f -f )(l-2«)du 
1 + u(l - u)*; 2 : 



so that it only remains to set z = 1 — 2u to derive (2.23). This concludes 
the proof of Lemma 4.3. □ 

4.3. The fourth order term. Our goal is now to provide an estimate on 
the fourth order term T 4 (A). We have estimated this term in (3.18), and 
we know that T 4 (A) is trace-class when A £ L 4 (M 3 ,IR 4 ). Here, we want to 
get an estimate involving only the norm of A in H^ iv { 



Lemma 4.4 (Estimate for the fourth order term). Let A = (A,V) G 

L 4 (M 3 ,M 4 ) n i^div( R3 ) and set 1 
universal constant K such that 



L 4 (M 3 ,M 4 ) n #div( R3 ) and set B ■= cullA and E '■= ~ W - There exists a 



| tr (tr C 4 T 4 (A)) | = | trT 4 (A)| < K (^ 1 -^) {\\B\\ L 2 + ||£|| L2 ) 4 . (4.19) 

V j=0 m jJ 

Proof. Arguing as in the proof of Proposition 3.1 (see the proof of Formu- 
las (3.31) and (3.34)), we decompose T 4 (A) as 

T 4 (A) = T 4)1 (A) - T 4 , 2 (A) + T 4 , 3 (A), (4.20) 

where 

T 4il {A) := ^- [ Cj u 2 du ^ \ 9(^2-5 \ z) , (4.21) 

2 W K ^ p 2 + m 2 +uj 2 \ p 2 + m 2 + uj 2 J K 

2 



i r v—-\ i / i 

27r JR f^ p 2 + m 2 + uj 2 \ p 2 + m 2 +uj 2 

x (Wi „ 2— 2 + ( Wi 2 W 2 -3— - ^ 



+ Wi = , V „ W 2 , ' = Wi 



p2 _|_ m 2 _|_ ^2 ^2 _|_ m 2 _|_ ^2 p 2 + m 2 + £J 2 _ 

(4.22) 
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and 

T 4 , 3 (A):=^ / Y^ Cj uj 2 doj 1 (wi 2 ■ \ t 2 ) . (4.23) 

2vr J R j-^ p 2 +mj + u 2 \ p 2 + mj + u 2 J 

Here, we have set for shortness, 

m := {p, A - Va} R3 + B ■ E - 2m, V/3, 

and 

W 2 :=|A| 2 + y 2 -2(a-^)y. 

Let us now explain our method to establish (4.19). When looking at Ta,u 
with k = 1,2,3, we are worried about several terms. First the function 
W2 does not decay too fast, it is only in L 3 (IR 3 ) if we only want to use 
the L 6 norm of A. Furthermore, it involves quantities which are not gauge 
invariant. Similarly, the term involving p in Wi involves non- gauge invariant 
quantities. On the other hand, the term involving B is in L 2 (M 3 ) and it is 
gauge invariant. The term involving V alone is also not gauge invariant but 
it has the matrix (3 which will help us, and it has no p. Since the result 
should be gauge invariant, these terms cannot be a problem. They should 
not contribute to the total (fourth order) energy. 

In order to see this, we use the following technique. In Formulas (4.21)- 
(4.23) , we commute all the operators involving p in order to place them either 
completely on the left or completely on the right. We have to commute the 
terms (p 2 + m 2 + uj 2 )~ 1 as well as the p appearing in W\. We think that it 
does not matter how many terms we put on the left and on the right. The 
main point is to have some functions of p on both sides (to get a trace-class 
operator under suitable assumptions on A). All the commutators obtained 
by these manipulations are better behaved and they will be estimated using 
the Kato-Seiler-Simon inequality (2.5), only in terms of ||A||^-i ^ R3 ^. 

In the end of the process, we will be left with a sum of terms of the form 

\P\ C f(x) \P\ d 

(p 2 + m 2 + uj 2 Y (p 2 + m 2 + u 2 ) b ' 

where f(x) is Wf or a product of W2 with some of the functions appearing 
in Wi, or only these functions. For instance, when we take the trace, the 
worst term involving only V is 

V*\[ I ,„ % ^ +3 f W 2dp 



(p 2 + m 2 + uj 2 ) 3 J K 3 (p 2 + m 2 + oj 2 ) a 



+ 



j 

4 



\p\ dp 



K s (p 2 +m 2 + u 2 } 



Here the integrals over p come respectively from 74,1, 74,2 and 74,3 and they 
behave exactly like (a; 2 + m 2 ) -5 / 2 . So we run into problems when we want 
to multiply by to 2 and then integrate with respect to uj. But this term 
cannot be a problem here because f^ 3 V 4 is not a gauge invariant quantity. 
This is where the Pauli-Villars scheme helps us. Not only these integrals 
will become well-defined, but also their sum will simply vanish because the 
regularization was precisely designed to preserve gauge invariance. 
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But before we explain all this in details, let us indicate how to handle the 
multiple commutators that we get when commuting the operators involving 
p. We start with T^i, for instance. Following the general strategy explained 
above, we write 



1 f 2 



+ 



3=0 
1 



(W : 



1 



p 2 + m 2 + uj 2 (p 2 + m 2 + oj 2 ) 2 



p 2 + m 2 + oj 2 



1 



p 2 + m 2 - + uj 2 



W 2 



3 3 

where, as usual, [S,T] := ST - TS. We notice that 



p 2 + m 2 + w 2 



W 2 , 



j,2 _|_ m 2 _|_ ^2 



-2- 2— [p 2 ,W 2 ] -j— 2 



p 2 + m 2 + uj 2 ' 



while 



b 2 ,w 2 ] = P [p,w 2 ] + [p,w 2 ] p = -i{p 1 vw 2 ] l 

Hence, we have 



1 r 2 

T 4,i(^) = ^- / Y. c ^ 2duj 

27T j=0 



1 



p 2 + m 2 + cj 



2 ^ 2 ^ (p 2 + to 2 + w 2 ) 2 



+ i 



Wo 



p 2 + m 2 + ui 2 p 2 + m 2 + Lxj 2 



{^VW 2 } R3 



(p 2 + m 2 + w 2 ) 2 



(4.24) 



where 



VW 2 = 2A • VA + W - 2V (a • VA) - 2 W (a ■ A) . 



We then argue as in the proof of Proposition 3.1. We use that W 2 G L 3 (M 3 ), 
with 



||W 2 || L3 <K||A|| 2 6 <tf||A||^ 



(R3)> 



and that VW 2 G L2 (R- 1 ), with 



||VW 2 || L 3 < K||A|| L 6||A||^ iv(ffi3) < K\\A\\%^ 

by the Sobolev inequality. By the Kato-Seiler-Simon inequality (2.5), we 
obtain for the term involving p ■ VW 2 , 

1 ... 1 1 



2 , 2 , -2 W2 — t 2— — P' VW 2 

_)_ m z _J_ ^z _|_ m t _|_ ^z 



(p 2 + m 2 - + a; 2 ) 2 
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< K 



p 2 + to 2 + w 2 



|VW 2 



w 2 
1 



6 3 



p 2 + m 2 + w 2 



|VW 2 



6 6 



(p 2 + m 2 + co 2 ) 2 
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< 



K 



(m 2 + uj 2 ) 



-A 1 1 si 



2^2 ll^lltf. 
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for some universal constant K. The argument is exactly the same for the 
term involving VW2 -p instead of p ■ VW2. Therefore, we obtain the bound 



/ 



2 

Em 

3=0 



1 



(p 2 + m 2 + w 2 ) 3 



1 



p 2 + m 2 + uj 2 



©I 



uj 2 duj 



<K 

V — ' m,- 



j=o 



I -All rrl 



4 



In particular, we have shown that the operator T 4i i(A) can be written in 
the form 



T 4 , 1 (A) = r 4 ,i(A) + 54,i(A), 



(4.25) 



with 



and 



l|r4,i(A)|| ei <^(X:M]|| A ||^ (ffi3) , (4.26) 



^4,i(A) /E« 



-J ^2 



+ rrij + or 



(W 2 ) 2 



(p 2 + m 2 + uj 2 ) 2 



uj 2 duj. 



(4.27) 

By the Kato-Seiler-Simon inequality, this term is trace-class when A € 
L 4 (M 3 ,]R 4 ) and conditions (2.12) are fulfilled. On the other hand, there 
is no evidence that the trace-class norm of 5 41 (A) can be bounded using 
only the norm ||VA|| i 2. Fortunately, this term will cancel with the other 
ones of the same type, as we will explain later. 

Our strategy to handle the operators T 4) 2(A) and T^(A) follows exactly 
the same lines. We first simplify the expressions of T 4 2 (A) and T 4 3 (A) by 
discarding the terms containing the operator B • "E. Concerning T 4i 3(A), we 
can compute 



1 



p 2 + m 2 + uj 2 



Wi 



1 



p 2 + m 2 + oj 2 



< 



K 



(m 2 + uj 2 ) 2 



\B\ 



L 2 



(||VA| 



3 4- 
L2 + 



m 



B ■ £ 

3 



1 



p 2 + m 2 + uj 2 
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3 



(m 2 + uj 



2\| 



|VF||1 2 + ||-B||| 2 ), 



so that 



27r -'K J^q P + m j + U V P + m ' 



with 



v .7=0 J 7 



(4.28) 
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We next commute, as above, the operator W\ — B • with the operator 
l/(p 2 + m 2 + oj 2 ) in order to establish that 

-!- / u 2 duy Cj - 9 \ ? f(Wi-5-E) - 9 \ 

27W R rr? V + mf + oj 2 \ ; p 2 + m 2 + w 2 / 

1 /• 2 
= r 4 6 3 (A) + — / ^ 2 ^^c,x 



, p z + + 



^2 j_ j_ /(1 2 

where 7^ 3 (A) also satisfies (4.28). Finally, we use that 

{p, A - a ■ V} R3 = 2p ■ (A - a ■ V) - ia ■ W, 

as well as the anti-commutation formulas for the matrices and /3, to 
obtain the formula 



1 f 2 ^ 

OJ 2 (lu> > Co — -jr k 7—r 

2vrA ^ 3 {p 2 + m 2 + u; 2 ) 4 



2 

/ W 2 liiJ 

({p, A - V a } RS - 2m,VpY p2 + ^ + uj2 = T 4 %(A) + 54,3(A), 



x 



with 7£ 3 (A) satisfying again (4.28), and 



2 

5 4)3 (A) :=- / ,^Vc )71 \ ((p 2 + m 2 ) 2 V 4 

3 

- 4(p 2 + m 2 )(m jy 9 + p • a)(p • A)V 3 /3 + 6(p 2 + m 2 ) • A)A l V 2 

i=i 

3 3 

- 4mj Yl PiPm(mj/3+p- a)(p • A)AiA m V 

1=1 m=l 
3 3 3 



1=1 m=l n=l / 1 



(4.29) 

The computation leading to this formula is tedious but elementary. In con- 
clusion, setting 74.3(A) := 743(A) +7^ 3 (A) + Tf 3 (A), we have established 
that 

74.3(A) = r 4 , 3 (A) + 5 4 , 3 (A) , (4.30) 
where 74.3(A) satisfies (4.28). Similarly, one can check that 

74.2(A) = 74.2(A) + 5 4 ,2(A), (4.31) 



with 



\T^A)\\ 6l <K(±m ||A||i L(ffi3) , (4.32) 
v .7=0 3 7 
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and 



S 4) 2 (A) 



- f Lo 2 duy"c jT ^ \ ^((p 2 + m 2 j )(3\A\ 2 + 

Wr ^ J (p 2 + m] + w 2 ) 3 V J V 



2(a • ^))y 2 - 2(p ■ a + m j p){p ■ ^)(3|A| 2 + 5V 2 )V 



+ 



3 £>(p ■ A)Ai (\A\ 2 + 5V 2 - 2(a • AJV)) 1 2 - 



(4.33) 



Notice here again that the Kato-Seiler-Simon inequality implies that £4,2 (A) 
and 54.3(A) are trace-class when A G L 4 (M 3 ,R 4 ) and conditions (2.12) are 
satisfied. Therefore, we always assume that A G L 4 (M 3 ,M 4 ) to make our 
calculations meaningful. 

The last step in the proof is to compute the traces of the singular operators 
5 4 ,i(A), 5 4 , 2 (A) and 5 4 , 3 (A) for A G L 4 (M 3 ,M 4 ). As announced before we 
claim that 



an identity which is enough to complete the proof of Lemma 4.4. To prove 
this we could make up an abstract argument based on gauge invariance. 
However we have to be careful with the fact that even if we can freely 
exchange the trace with the integration over uj, these only make sense after 
we have taken the sum over the coefficients Cj. The order matters and 
this complicates the mathematical analysis. Instead, we calculate the sum 
explicitly and verify that it is equal to 0. 

A simple computation in Fourier space shows that 



trS 4 ,i(A) - tr 5 4l2 (A) + tr 54.3(A) = 0, 



(4.34) 




(4.35) 



x 



it 3 



/ (\A\ 4 + 6\A\ 2 V 2 + V 4 ) 



Similarly, one can check that 





3 3 



An integration by parts shows that 
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and we obtain 
tr<S 4 , 2 (A) 

„2 , ™2 



= -^f dp [ ^fe Cj - / + 2 m ' 2U ) / (MV + v*) 

+ /" d p /" u^fy-cj— \ [ (\A\* + 5\A\ 2 V 2 ). 



j=o v ^ J 

(4.36) 



Similar computations lead to the expression 

32 / /■ /• 2 (p 2 +m 2 ) 2 rfp^d^ f 

+— ^— r / / y c d ^ 2d " ri i 3(p2 + m2) ^ /" i i|2 T - 2 
vr(2vr)iURA3^ Cj (^ + m 2 + w2) 3^ + p2 + ro 2 + w2 ;; y ffi3 1 1 

2 ( f [ dpu 2 du \ I" 4 

(4.37) 

In view of (4.35) and (4.36), we obtain 
trS 4 i(A) -tr<S 42 (A) + trS 43 (A) = — — T ( / V 1 ) / oj 2 doj / dpx 



(27 

x Tc ( 1 - 12 p2+ ^ 2 + i 6 (P 2 + m ?) 2 

^ 5 V (p 2 + m 2 + w 2 ) 3 (p 2 + m 2 + w 2 ) 4 (p 2 + m 2 + w 2 ) 5 



A direct computation then shows that 

1 12 16 \ 2 , 

This is enough to deduce (4.34), and complete the proof of Lemma 4.4. □ 

4.4. Regularity of the sixth order term. In this section, we come back 
to the sixth order term studied in the proof of Proposition 3.1. The sixth 
order term is defined as 

K 6 {A) = -L f tr (R' 6 (lo, A) + R' 6 (-u, A)) doo, (4.38) 
4tt J r 

where 

BL{u,A) :=y Cj l - Uoc-A-V) Y. 

We have shown that it is trace-class when A G (M. 3 ). We can indeed 
write estimate (3.14) as 

2 

H^^lk ^ ^ l I W w > A)\\ 6i du,<K(J^W) \\A\\% L{R3y (4.39) 

3 — ^ 
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Here we want to prove that TZq is actually smooth, under suitable assump- 
tions on A. We first establish the continuity of TZq through 

Lemma 4.5 (Continuity of the sixth order term). The functional TZq is 
locally 9-Hdlder continuous on the space ijJ[ iv (IR 3 ) for any < < 1. 

Proof. We consider the difference R' 6 (uj,A) — R' 6 (u,A') for four-potentials 
A and A' in a given ball of ij^ iv (IR 3 ), which we write as 



R' 6 (uj,A)-R' 6 (u,A') 

2 

3=0 
5 

' l),„ 4 i 

k=0 1 ' 



1 



D mj ,A + iu D m . iA , + iuj 



a 



■ A — V) 



(( a -A'-V)—± Y: 



(a • (A - A') - V + V') 



D mi ,o + iu 



a 



, 1 \ 5 - k \ 
■A-V)— 

D mj ,o + iuJ J 

(4.40) 



The five terms in the sum over the index k can be estimated similarly as in 
the proof of (3.14). Their ©i-norms are bounded by a universal constant 
K times 



j=o (m 2 j +uj 2 )2 



+ \\A 



/i|5 



I A An,, 



^iv( R3 )' 



If we follow the same proof for the first term, we need an estimate on the 
operator norm 

1 



1 



- D mj , a + D mj iA > + iw' 
On one hand, we remark that \\(D m a + w) _1 || < 1/oj, so that 

1 



iuj 



( l - 



D mn , A + iu D m . Ai+iu* 



< 2. 



On the other hand, we can use the resolvent formula to write 
1 1 



D mjj A + iu D m . :A >+iu 



D m ,, A + iu 



(a -(A- A') + V'-V) 



1 



D m ,, A' + W 



(4.41) 



(4.42) 



For small A, or small A', we have no problem in estimating this term using 
that the spectrum of D mj ,A stays away from by Lemma 2.1, and that 
(D mj ,A + iw) _1 (D m ^o + ioj) is bounded uniformly. The argument is essen- 
tially the same in the general case. We decompose the expression in the 
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right-hand side of (4.42) as 
1 



iuj 



D mi ,A + iuJ 



[a - (A — A') + V- V) 1 ) 



= (-=. 1 ——{D m A + in)) x (- 1 ——(D m + iv)) 



( a .(A-A') + V'-V)) 



D m ,, A' + iu' 



(4.43) 



for some positive number [i. We check that 



(D mj ,A + ilj) 



< 



An, ,A 



D mj>A + iuj 



+ 



i/j, 



D mj>A + iu 



< 1 + 



Setting a := 4i^ 2 ||A||^ >1 , , we also remark that 



\U\ 

(4.44) 



Miv(R 3 )' 



Dm^A + ifJ- 



1 



(D mj ,o + = (1 + (V - a ■ A 



oo -. 

<y II— — — (v-a-A) 



(4.45) 



K r ' 



— II 411™ < 2 

n ll^-ll III (TB>3\ 



„=oK 2 + / i2 ) : Hdiv 

Recalling that {ioj)\\{D m ^ AI < 1, we infer from (4.42), (4.43), (4.44) 
and (4.45) that 



(iu) (— : — — ) 

K 



< 



{m) + ii 2 ) 4 



if 1 + 



(4.46) 



<K 



1 1 I, . I, 

fPPT + I, ,1 H^Hh 1 . (R3) 



\A-A'\\ Hl 



In this bound, we can replace A by A' by symmetry. Recall that we are in 
a given ball in i^ iv (M 3 ), so that ||A||^i ^ R3 ^ is bounded by some constant. 

Collecting estimates (4.41) and (4.46), we have shown that 
/ \\B! & {u,A)-B! & {u,A')\\ &1 cko 



I 4 II 5 -I- II 4'll 5 



+ K ( i2 j^f j (H A H^ iv (R3) + II A/ |l^ iv (R3))^K-> A, A'), 



A A 'll^ iv (R 3 ) (4.47) 
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/>oo 

I( mj ,A,A'):= - 
Jo (1 



doj 



+ W 2 )2 



x mm 



I 2 ' (,/m- + m lW l|A|l ^iv(^ 3 ))ll A A/ H^ iv (^ 3 )}- 



Assuming that \\A - A'H^i (|3) < rrij/(^/fnj + ||A||^-i , R3 s) for any j = 



0, 1, 2, we can estimate the integral l(mj,A, A') as 



Pto,A,A')\<-jL:(l 



with 



(1 + w 

Jo (1 + 0J 2 )2 Jt 



\A-A>\\ m 



\A-A'\\ Hl 

dr 



dco 



3 * 

>t Uj{l+0J 2 )2 



It remains to observe that 



J{t)<Kt{l + \\ogt\), 

and to combine with (4.38) and (4.47), to conclude that the functional H§ 
is locally <9-H61der for any < 9 < 1. □ 

We next turn to the differentiability of TZq. 

Lemma 4.6 (Regularity of the sixth order term). The functional TZq is of 
class C°° on the open subset H of H^ iv (]R 3 ) containing all the four-potentials 
A such that ^ o-(D mj ,A) for each j = 0,1,2. Moreover, there exists a 
universal constant K such that 



|d 2 ^ 6 (A)||<K^M 1 + ^||A||^ 

3=0 j 



(4.48) 



where 



L A --max {\\(D mj:A + iLo) 1 (D mjt0 + u€R, j = 0, 1, 2} < oo. 

Proof. The proof relies on elements in the proof of Lemma 4.5. When is not 
an eigenvalue of D mjt A for each j = 0, 1, 2, we can deduce from Lemma 2.1 
the existence of a positive constant Ka such that 



1 



D m ,, A' + iu 



< min 



(4.49) 



for any A' 6 i^ iv (IR 3 ), with \\A' — A\\^i , r3 n small enough. As a conse- 
quence, we can replace estimate (4.44) by the inequality 

1 



D mj ,A + iu 



< 1 + K Af i. 
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Since 



1 

D m -A + iu 



it follows that 
1 



(D mj fi + iu) 



< 



D nijtA + iuj 



{D mj! o + i[i) 



+ \u — n\ min 



D m ,A + iu 



(D mjfi + iu) 



< (1 + K^) 



Dm^A + ifJ- 



+ | a; — /x| min j-fCl, j — j- j- 

(4.50) 



Following the lines of the proof of (4.45), we deduce that the quantity in 
the right-hand side of (4.50) is bounded independently on u by a positive 
constant La, depending only on the four-potential A and the mass rrij. 
Actually, we can claim, up to a possible larger choice of La, that 



D m .A> + iu 



(D mjfi + iu) 



< L A , 



for any u G M, j = 0, 1, 2, and A' G i^ iv (M 3 ), with || A' - A|| (R3) small 
enough. 

As a result, we can upgrade (4.46) into 



j[ I 

^ D mj ,a + iu D m . a 1 +iuJ 

Similarly, we can compute 



< KL * r \\A-A'\ 

(m 2 + u 2 ) 4 



IU 



D mj>A + iu 



(a -(A- A') + V-V) 



1 



D m ,, A + iu 



KL A I, ,m 



(4.51) 



(m^ + u 2 )i 



At this stage, we can iterate the resolvent expansion in (4.42) to obtain 
1 1 



D m . A + iu D m . A '+iu 



(a • (A - A') + V-V) ' 



An^A + iu 



D mi) A + 



+ (~R (a • - A') + y' - F)f 



mj,A' + 



Inserting this identity into (4.40) , we can write 

R' 6 (u, A) - R' e (u, A') = d A R' 6 (u, A)(A - A') + r' 6 (u, A, A') . (4.52) 
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Here, d^i^^w, A) refers to the continuous linear mapping fromi^QR 3 ) to 
6i(M 3 ,R 4 ) given by 



d A R' 6 (u, A)(t),o) 

= I>(V^ (a • a — t>) — ((a-A-V)— — V 



icj / . In 5-fc 



+ E n ■ (»^^) n • ( (a ■ a — e) — — 



fc 

+ ioj 



(4.53) 

In view of (4.51), and again the computations in the proof of estimate (3.14), 
the operator norm of d^-Rg^, A) is indeed bounded by 

||d A i?6KA)|| 

2 i I 

- K § k +- 2 ) f l|i4fe - (R3) ( 1 + K +^)i l|A|l ^v(« 3 ))- 

(4.54) 

Similarly, the remainder r 6 (w, A, A') in (4.52) may be estimated as 
r 6 (u,,A,A') 

6i 

2 ^ | 

< K|| A - A'\\% Lm g ( 2 + J w2) | (II A H^ iv (R3) + H^'fe^CRS)) X 
^2 

X ^ + (m 2 +l 2 )^ l|A|l ^(^) + l|A/|l ^ iv (K 3 )))" 

(4.55) 

Collecting (4.38), (4.52), (4.54) and (4.55) is enough to establish the con- 
tinuous differentiability of the function TZq on a neighborhood of A, with a 
differential given by 

dft 6 (A)(o,a) = / tr fd^e^, A)(o,a) + d A i? 6 (-w, A)(o, a)) dw. 

4vr J M V / 

(4.56) 

Finally, we can extend the previous arguments for the continuous differ- 
entiability of TZq to the proof that it is actually of class C°°. In particular, 
we can check that the norm of the quadratic form d\R' 6 is bounded by 



d^i^A) 



L 



2 



— ^ ^ ^ 3 1 1 -^4. 1 1 jjl ^]J^3 ^ ( ^ ill 1 1 i 1 / ^ 

(m 2 + w 2 )5 divl - J n (m 2 + uj 2 )? div 
Estimate (4.48) follows integrating with respect to w. □ 
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When ||A||^i ^ R3 ^ is small enough, we can prove that the constant La 
does not depend on A. 

Corollary 4.1 (Estimate in a neighborhood of zero). Assume that Cj and 
rrij satisfy (2.15). There exists a universal constant rj such that, given any 
A 6 i7j iv (M 3 ) with \\A\\^i ^ ffi3 ^ < rj^/mo, the functional TZq is of class C°° 
on the ball 

B A {vVrno) = {At Hl iv (R 3 ) : \\A\\^ {R3) < n^/T^} , 
and satisfies the estimate 



|d 2 ^ 6 (A)||<;d ]T 



n m i 



(4.57) 



Proof. When A is small enough, the spectrum of D mjt A does not contain 
by Lemma 2.1. Moreover, when 



I A 1 1 Tjl 



for r] small enough, we can infer from (4.45) that 

L A < 2. 



Inserting in (4.48), and using the inequality ||A||^i ^ ffi3 ^ < r/^/mJ, gives 
estimate (4.57). □ 



5. Proof of Theorem 2.1 

With the results of the previous section at hand, the proof of Theorem 2.1 
is only a few lines. As a matter of fact, given any A £ L 1 (M 3 , M 4 )ni/j iv (lR 3 ), 
we have shown that the functional J r pv(A) is well-defined by the expression 

T PY (A)=T 2 {F) + K{A), (5.1) 

where 

F 2 {F) :=tr(tr c4 T 2 (A)), 

and 

K(A) := tr ( tr C 4 T 4 (A)) +tr(tr C 4T^(A)), 

are defined in (4.1). By Lemma 4.2, the function Ti is given by (2.20) and 
it is quadratic with respect to F. Since M is bounded, we deduce that F2 
is smooth on L 2 (M 3 ,IR 6 ). On the other hand, the function A i-> F^(A) := 
tr(tr C 4 T 4 (A)) is quartic and satisfies (4.19). Hence, it is a smooth func- 
tion on i5r^ iv (M 3 ). We have proved separately in Lemma 4.5 above that 
TZq (A) = tr(tr C 4 Tg(A)) is an Holder continuous function on i^ iv (IR 3 ), which 
satisfies (3.14). We deduce from all this that J-~py has a unique continuous 
extension to i^ iv (M 3 ), which is given by (5.1), and that 1Z satisfies esti- 
mate (2.19). The properties of M can be found in Lemma 4.3. □ 
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6. Proof of Theorem 2.2 

In view of the results in Sections 3 and 4, the functional Tpy is smooth 
on the open subset % of ijJ[ iv (R 3 ) containing all the four-potentials A such 
that ^ cr(D mj) A) for each j = 0, 1, 2. Indeed, the function A i— >■ T^{A) := 
tr(tr^4 T±{A)) is quartic and satisfies (4.19). Hence, it is of class C°° on 
i5r^ iv (]R 3 ). Similarly, in view of Lemmas 4.2 and 4.3, the quadratic map Ti 
is smooth on L 2 (M 3 ,IR 6 ). On the other hand, we have shown in Section 4.4 
that 7^6 is smooth when is not an eigenvalue of D mji A for j = 0, 1, 2. We 
deduce that Fp\ is smooth on the set %. 

In order to complete the proof of Theorem 2.2, it remains to identify 
dTpv(A). As mentioned in Formulas (2.24) and (2.25), this differential is 
related to the operator 

2 

QA ■= C i -"-(-00,0) {Prrij,A) ■ 

i=o 

Concerning the properties of the operator Q a , we can establish the following 

Lemma 6.1 (Properties of pa and ja)~ Assume that Cj and rtij satisfy 
conditions (2.15). 

(i) Let A £ -ff^ iv (M 3 ) be a four-potential such that is not an eigenvalue of 
D mj ,A for j = 0, 1, 2. Then the operators tr C 4 Qa and tr C 4 olQa are locally 
trace- class on L 2 (IR 3 ,IR 4 ). More precisely, given any function \ £ L^°(M 3 ) 
(that is, bounded with compact support), the maps 

AeH^ x(tr C 4 Q A )x G ©i 

and 

A^U^f x(tr C 4 c*Q A )x G ©i 
are continuous from H to &\. In particular, the density pa and the current 
j A, given by 

Pa{x) := [tr C 4 Qa] (x,x) and j A (x) := [tr C 4 a Q A ] (x, x), 

are well-defined and locally integrable on M 3 . Moreover, the maps A i— >■ pa X 2 
and A i-> j'a X 2 are continuous from H to L 1 (M 3 ). Finally, for A = 0, we 
have 

p = and j = 0. 

(m) // moreover A G L 1 (M 3 ,IR 4 ), i/ten, i/te operators tr C 4 (Qa — Qo) ffln ^ 
tr C 4 q((5a — <3o) « r e trace-class on L 2 (M 3 ,R 4 ), and the density pa and the 
current 3a are in L x ( 



Proof. We split the proof into three steps. First, we consider the special 
case A = 0. 

The operators tr C 4 Qq and tr C 4 olQq are locally trace-class. Using that 
Ylj=o c j = 0) we can write 

1\ 1 \r-^ D m .n 



Qo = J2cj (l(-oo,0)(Ary,o) - g) = "^ tr C 4 X^ 



j=0 "' j=0 ' mj '° 
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As a consequence, we obtain 

tr C 4 Qq = 0. 

In particular, the density po is well-defined and it identically vanishes on M 3 . 
Similarly, we have 

2 



tr C 4 aQ = " 



Due to conditions (2.12), the latter function behaves like 

2 „ 2 

E P 3 ^ 4 g 

Cj ( 2 i 2^ ~ a Z^ c J m i |„|5' 

i=0 (p 2 + mj) 2 » i=Q |P| 

as |p| — > oo, hence it is in L 1 (M 3 ). By the Kato-Seiler-Simon inequality (2.5), 
we deduce that the operator tr C 4 axQoX is trace-class for any x £ L 2 (M. 3 ). 
Hence, tr C 4 aQo is locally trace-class. In particular, the current jo is well- 
defined and locally integrable on M 3 . Moreover, we can compute 



/ Jo X 2 = tr ( tr C 4 axQox) 

JR3 



7^3 / / J2c j —^-— T \x(q-p)\ 2 dpdq, 



which shows that 



1 f 2 

jo = - -r^T / / Cj r dp = 

4W R 3^ J (p 2 + m^ 

by rotational symmetry. 

We next consider the general case. 

The operators tr C 4 Qa and tr C 4 olQa are (locally) trace-class. From 
the previous discussion, we conclude that it is sufficient to prove that the 
operators ^^{Qa — Qo) an d tr C 4 ol{Qa — Qo) are locally trace-class. The 
corresponding charge and current densities will be the same as that of Qa- 
Concerning the (local) trace-class nature of the operator Qa — Qo, we 
follow the proof of Proposition 3.1. Our starting point is the integral formula 

2 f xoj 2 1 f ( iu ioo \ 

When T is a self-adjoint operator on L 2 (R 3 ,R 4 ) with domain D(T), we 
deduce that the sign of T is given by 

signT = — / - du>, 6.2 

the integral in the right-hand side of (6.2) being convergent as an operator 
from D(T) to L 2 (M 3 ,C 4 ). 
In particular, the operator 

! 2 

Qa - Qo = 5Z C J ( si S n An.,-,A - signD m3i o), 
i=o 
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is given by the expression 
i r 2 

1 / \ - / UxJ ILU 



+ 77^ '■ \9 



{D mj , A -iu) 2 (D mjt0 - iuj) 2 , 

on i? 1 (R 3 ,C 4 ). In order to establish statements (n) and (ra) of Lemma 6.1, 
we will prove that 

2 



(An, ,A + iu) 2 {D m . >0 + iw) 

da; < oo, 



lUJ lUJ 

+ 7~R 7-v> X 



{Dm 3 ,A-iu) 2 (D m0 -iuj) 



Si 

(6.4) 



for any of the matrices m = I4, a±, 0:2, 0:3, and either when A G L^M 3 ,!^ 4 )!"! 
^(IR 3 ,^ 4 ) and x = 1, or when A G #div( M3 ) and X G L2°(M 3 ,]R). In the 
different cases, the C 4 ~traces of the operators m(QA — Qo)j respectively 
wx(Qa — Qo)Xj will define trace-class operators on L 2 (M 3 ,C 4 ). Then the 
operators tr C 4 vciQa will be locally trace-class and the density pa and the 
current ja will be well-defined and locally integrable on M 3 . Moreover, they 
will be integrable on M 3 for any A € L^IR 3 ,^ 4 ) n tf^R^R 4 ). 
In order to prove (6.4), we use the expansion 

ioj ioj ioj iuj 

+ 



(D mj ,A + iuj) 2 (D mjt0 + iu) 2 {D mj! A ~ iu) 2 (D mjfi - °> • 



:= (Qn(u, A) + Q n (-co, A)) + Q' 6 (u, A) + Q' e (-co, A) 
- Q' 7 (uj,A) - Q' 7 (-u,A), 



VjJ ) 

(6.5) 



n=l 



with 



2 In 

QJu, A) := (n + 1) V Cj — ^ f (a ■ A - V) Y , 

Q' 6 (u, A) := 7 V cj — iUJ . - 2 ( (a ■ A - V) - 1 , 



3=0 

and 

2 

Q' 7 (w,A) :=6y c?7 -f(a-A-F) ) . 

We next estimate the terms related to the operators Q n (uj,A), Q'q(uj,A) 
and Q' 7 (u!,A), as we have previously done for the operators R n (co,A) and 
R'q(uj, A) in Section 3. 

Concerning Q' 6 (w, A) and Q 7 (uj, A), we recall that is not an eigenvalue 
of D m .A for each j = 0, 1, 2. Hence, there exists a positive constant K such 
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that 



D mj t A + iw 



<K, 



(6.6) 



for all u 6 R and j = 0, 1, 2. Following the proof of (3.14), we deduce that 
(\\Q' Q {u3,A)\\ &1 + \\Q' 7 {u,A)\\ &1 )dw 

1 



/ 



n m 7 



21 HV4ll 6 o 

2 H Vyl llL 2 



1 + 



/rrii 



IVAI 



L 2 



As a consequence, the integrals 

Q' 6 (A) := i- jf (Q^( W , A) + £%(- w , ^)) d", 



and 



Q' 7 (A) :=-L / (Q 7 ( W ,A) + Q 7 (-^,A))da;, 



when A G if] iv 



The related 



define trace-class operators on l?{ 
densities p' e (A) and p' 7 (A), and currents j' e (A) and j 7 (A), are well-defined 
and integrable on R 3 . Moreover, in view of (6.6), we can repeat the argu- 
ments in the proof of Lemma 4.6 in order to establish the smoothness of the 
maps A i-> Q'q(A) and A i-> Q' 7 (A) from H onto 6i. 

For 3 < n < 5, the operators Q n (^>, A) satisfy the estimates 



and 



||Q n (w,i4)|| 6i dw<if n ||A| 



/ ||xQn(w,i4)x|| 6l duj < K n \\A\\ n L6 \\x\\ 2 LT M-, 



(6.7) 



(6.8) 



for any function x £ L£°(R ). Here, K n refers to a positive constant de- 
pending only on the coefficients Cj and the masses rrij. For n = 4 and n = 5, 
we can indeed use the Kato-Seiler-Simon inequality (2.5) to write 

2 



\Q n {u,A)\\ 6i KK^c^Af^ 



|w| dp 



i=0 (p 2 + m| + w 2 ) 2 

Integrating with respect to w, we obtain inequality (6.7) with 

n n-3 
3=0 m 3 

For n = 3, we rely on the identity Co + c\ + C2 = to write 
Q 3 ( W ,A) = 4y c / V-^ ( — — (a-A-V)) k x 



3=0 



( 



+ 



A-V) 



3-k 



D mjfi + iu D mofi + iuj 
iuj iuj 



a 



D mj ,o + iu< 
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Using inequality (3.16), we deduce that 

|3 



Qs(u,A)\\ < K[ J^lcjKmj - m )) ^ . ||A|| L ,. 

V „_ n / (mi +oj 2 )2 



which provides estimate (6.7) with 



,"=0 



Inequalities (6.8) follow similarly. Applying the Sobolev inequality (3.12) 
to (6.8), we deduce that the integrals 

Q n (A) :=^- I (Q n (u,A) + Q n (-co,A))duj, 

define locally trace-class operators on L 2 (IR 3 ,]R 4 ) for 3 < n < 5, as soon as 
A G i?^ iv (lR 3 ). The related densities p n (A) and currents j n (A) are well- 
defined and locally integrable on R 3 . When A is moreover in L"(1R 3 ), in- 
equality (6.7) guarantees that the operators Q n (A) are trace-class, while the 
functions p n (A) and j n (A) are integrable on M 3 . The continuity in these 
spaces follows from multi-linearity. 

For n = 1, we refine our estimates using the cancellations provided by con- 
ditions (2.12). Following the lines of the analysis of the operator R±(uj, A), 
we start by writing 



Qi(w, A) + Qi(-w, A) = Q 1A (u, A) - Q lj2 (w, A), (6.9) 

where 

~q K-^rrij ,0 

(6.10) 

and 



Q lA u,A) ^ {a • A - ^, } R3 ^_ : 



Q 1>2 ( W> A) := 4 £ c, a; 2 (a ■ A - V) ■ 



As for the operator Qi,2(w,A), we combine conditions (2.12) with identi- 
ties (3.26) to estimate 

||Q 1 , 2 (a;,A)|| 6i <K(^|c J |(m 2 -m 2 ) 2 ) f \\A\\ LX . (6.11) 

V „--n / (mi +u z 2 



i=o 



In order to estimate the operator A), we eliminate the odd powers 

of the masses mj in the numerator of the right-hand side of (6.10) by taking 
the C 4 -trace. Recall that 

{a- A- V,D mjfi } R3 = {p,A- Va} R3 + B ■ E - 2m,V73. 

Since 

tr C 4^n^] =0, (6.12) 
^ k=i ' 
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when d is odd, we obtain 



tr C 4 (mQi,i(w,j4)) :=8^ Cj a; 2 ——- ( - 2m 2 Vtr C 4 (m) 

+ tr c4 (m(a • p) ({p, A - Ka} R3 + 5 • -— L 



+ 0J 2 

for any of the matrices m = c*i, ct; 2 , 0:3. On the other hand, we can 
compute 



1 



1 



+ 



2 2 
m,Q — m- 



(p 2 + m 2 + lo 2 ) 2 (p 2 + rriQ + uj 2 ) 2 (p 2 + m 2 + LJ 2 )(p 2 + rriQ + uj 2 ) 2 

2 2 
777.Q — m- 



(p 2 + m 2 + oj 2 ) 2 (p 2 + ijiq + w 2 ) 



(6.13) 



as well as 



2{m\ - m 2 ) 



(p 2 + m 2 - + uj 2 ) 2 (p 2 +wig + w 2 ) 2 (p 2 + m 2 , + ''' 2N > 3 



+ 



2(?tiq — to 2 ) 2 (ffiQ - m 2 ) 2 



(p 2 + m 2 + uj 2 )(p 2 + to 2 , + w 2 ) 3 (p 2 + to 2 + w 2 ) 2 (p 2 + m 2 , + a; 2 ) 2 

Combining again with conditions (2.12) and identities (3.26), we obtain the 
estimate 

2 ^2 

tr C 4 (mQi,i(w,j4)) < | Cj | (m 2 - m 2 ) — -yX 

- , „ (m5 + o; 2 )2 



j=0 

x 



(m 2 ||y|| L1 + (m 2 -m2)||A|| L1 ). 



In view of (6.9) and (6.11), we have 

/ tr C 4 (m (Qi(u, A) + Qi(-u, A) 



61 



duj 



<Kj2\cj\H-m 2 ) 

j=0 

Similarly, we can check that 



mj — m 



- \\v\\ , + J 

'"-o 







4 

2 ii^IIl 1 



(6.14) 



/ tr C 4 (mx(Qi(w,A) + Q 1 (-u,A))x) 



61 



da; 



i=o 



< Kj2\ci\ K - m 2 ) ( -I ||y|| L6 + ||A 



lL6yi|X|| L f- 



(6.15) 



For n = 2, the analysis is identical. We compute 

Q 2 (w, A) + Q 2 (-w, A) = Q 2jl (w, A) - Q 2)2 (w, A), (6.16) 
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where 

2 

Q 2 ,i(u,A) := 12J2 c J uj2x 

3=0 



and 

2 



+ 7^ 2-{ a - A - V ) r,2 Dmj '° 2 i a - A - V )-F^ 2 



+ L 2 1 2 (a • ^4 - V)) \ 2 Dm »° 2 ) 



In order to estimate the operators Q 2) i(o;, A) and Q2,2^, A), we again take 
the C 4 -trace. For m = I4 or m = a.j,, we derive from (6.12) that 



2 

2 



tr C 4 (niQ2,i(w,A)) =12^c 



i=o 



-2m 2 m— r ^ ^V-^ Act-A-V) ° 



a • p jr 1 , _ A i 1 



3 iPLfi+^Y d 2 +^ >D 2 + 



p 2 + m 2 1 1 \ 
+ — o a ■ A-V.a - x? V — ^ -(a • A - V)—^ . 

while 



tr C 4 (mft,2(w,A)) = 6^Cja; 2 tr C 4 (2m 



X 



((a-A-V)—*-^ »Y +m^ r -^ -(a- A-V) 

(a- A- V)— 5—^ 5 + m(— 5— ^ 7 (a ■ ^4 — V")) 



a • p 



2 a • p 

Dl^ + u 2 
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Invoking conditions (2.12), as well as identities (3.26) and (6.13), we deduce 
that 



tr C 4 (mQ 2 ,i(u>,A)) <^V|cj| ( (m 2 - m 2 ) \\A 

J=0 



IL 2 



(m,Q + oj 2 ) 2 

2 



+ m ) \\ A \\ L - (ii a \\v + imw (m2 ^ 2)2 ) > 



and 



tr C 4 (mQ 2 ,2(w, A)) < k( \ Cj \ (m 2 - m 2 ,) 
1 S=o 



lL2 (m 2 +w 2 ) 2 ' 



It follows that 



/ tr C 4 ( m(Q 2 (uj, A) + Q 2 (-uj, A) 

JR V 



6i 



Similarly, we have 



2 2 

m- — rriQ 
m 



\ A \\v+ m i\\ A \\v(\\ A \\l? + \\ V \\v) ■ 



tr C 4 (mx(Q2(w,A) + Q 2 (-w,A))x) 



©i 



2 / 2 2 

771- m^ 1 1 2 ii ii /ii ii n n \ \ 1 1 

— II^Le +mj \\A\\ L6 {\\A\\ L6 + \\V\\ L6 ) \\x\ lL3 . 



3=0 



m 



(6.17) 



In view of (6.14) and (6.15), we conclude that the integrals 

tr C 4 (m Q n ) := — tr C 4 (m (Q„(uj, A) + Q n (-u, A)) ) duj, 

also define local trace-class operators on L 2 (M 3 ,R 4 ) for n = 1,2, as soon 
as A £ H^ iv (R 3 ). The operators are trace-class when A is in L"(IR 3 ). 
Concerning the related densities p n (A) and currents j n (A), they are well- 
defined and locally integrable on M 3 for A G i^ iv (M 3 ), and integrable on M 3 
for A € L n (IR 3 ). Their continuity follows again by multi- linearity. 

At this stage, it remains to recall Formulas (6.3) and (6.5) to complete 
the proof of Lemma 6.1. □ 



We are now in position to complete the proof of Theorem 2.2. 

End of the proof of Theorem 2.2. We have shown that the functional Jpy 
is smooth on the open subset % of four-potentials A such that is not 
an eigenvalue of D mjt A for each j = 0,1,2. In particular, the differential 
dJ r pv(A) is a bounded form on i5fj iv (R 3 ). By duality, it can be identified 
with a couple of functions (p*, j*) in the Coulomb space C defined in (2.26). 
Our task reduces to verify that = pa and = —ja- 
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We first restrict our attention to four-potentials A which are moreover 
integrable on M 3 . In this case, the functional .Fpv(A) is given by For- 
mula (2.17), which may be written in view of (4.1) as 

5 

T PY (A) = J2MA) + U 6 (A), 

n=l 

where we recall that 

Fn(A) := ^^tr (tr C 4 (R n (uj, A) + R n (-uj, A))^j du, 



and 



TZq(A) := i- J tr ( tr C 4 (R' 6 (u, A) + R' 6 (-cu, A)) dw. 



We have computed the differential of &1Z§(A) in (4.56). On the other hand, 
the functionals T n are n-linear with respect to A, so that their differentials 
are given by 

cLF n (A)(o,o) 

= J tr(tr C 4 (d A R n (u, A)(e,a) + d A R n (-u, A)(t),a))) dw, 

with 

n-l 



(U-RJw, A)(o, a) = Vc 7 - — — V ((a ■ A - V)— — )*x 



n—l—k 

X 



(a • o - t»)— ((a ■ A - V)— ) 



(6.18) 



for any (t), a) G L 1 (E 3 , R 4 ) n #] iv (K 3 ). It follows that the differential 
d.Fpv(A) is equal to 

dFp V (A)(t),a) = — / S(w,i4)(t>,o)dw, 
4vr J R 



with 



/ 5 

S(w,i4)(t>,a) :=tr I J^tr (tr C 4 (d A i? n (w, A)(o, a) + d A R n (-uj, A)(tt, a))) 

+ tr (tr C 4 (dA# 6 (w, A)(D,a) + d A i?e(-w, A)(o,a)))V 
At this stage, we make use of Formulas (4.53) and (6.18) to check that 

E{ U ,A)(p,a) = tr (tr c4 ( g Cj ^ ™ + (a ■ a - *))). (6.19) 

Indeed, we have established in the course of Lemma 4.6 that each term in 
the decomposition of d A R' 6 (co, A)(t), a) which is provided by Formula (4.53) 
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is trace-class. As a consequence, we can write 
tr(tr C 4 d A R' 6 (co, A)(o,a)) 

= Vq tr — — — (a ■ a - t>)— ((a ■ A - V) — V 

4-i J D mi A+iuj y ' D m . A + icu\ K 'D m . + tu)J 



2 5 



+ 5> Vtr— -™ f(«.^-y)_J ) 5 ~\ 



j=0 fc=0 



x (a-o-t))- — ((a-A-V)- — Y . 

D mj ,o D mj:0 + iuJ 

An advantage of this further decomposition is that we are allowed to com- 
mute the products in the right-hand side, so as to obtain 

tr(tr C 4 d A R' 6 (u, A)(o,a)) 

2 

= Vc,tr— — — ((a-A-V)— ) — (a ■ a - t>) 

J^o D mj)A + iuV >D mjfi + iLj) D m]A + iu y > 

+ Tc- Vtr l - ((a A V) 1 Y 1 x 

~" 3 fr D mjfi + iujV ' D mj:0 + iujJ D mjtA + iio 

x ((a-A-V)- — Y~ k (a- a-t>). 

This follows from the property that the operator (ico)(D mjjA + iu)~ L is 
bounded, while the operators (a • A — V)(D mjj o + and (a ■ a — 

b)(D mj: o + iuj)^ 1 belong to suitable Schatten spaces. Using the resolvent 
expansion (3.6), we are led to 

tr (tr C 4dA#6(w, A)(D,o)j = tr ( tr C 4 ^ Cj ( — ^ . . 2 (a ■ a - t>) 

-VV-i ( {a . A -V)—^ Y - ™ x 

fclo t^O m ^'° + ^ D m 3 fi +IUJ D mjfi + tUJ 

(6.20) 

Similarly, we can deduce from (6.18) that 
tr (tr C 4 d A i? n (w, A)(r-,a)) 

(2 n— 1 ^ ^ n _~i_f t 
tr c4 Vc 7 V— C( a .A-V)— ^ x 
^ fro + w V + ^ 

x 7^ ((a-A-V)— )*(<*■ a- t>)Y 

Formula (6.19) follows combining with (6.20). 
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As a conclusion, we have derived the following expression of dJ r pv(A), 
i r / / 2 



dJp V (A)(t),a) = -L jTtr(tr C 4 ^ 



^ (a ■ a — t)) doo^j . 



(D mj ,A ~ iu) 2 
In view of (6.1) and Lemma 6.1, we deduce that 

dJ"p V (A)(o,a) = tr (tr C 4 (q a {*> - a ■ a))) = J ^ {pa^-Ja ■ a), 



so that = pa and j* = j^, when A G L 1 (R 3 ,R 4 ) n -ff^ivV 

In the general case where A is only in ij^ iv (IR 3 ), we can construct a 
sequence of maps (A Tt ) ngN in L X (IR 3 , ]R 4 )ni?] iv (IR 3 ), for which ^ a(D mjtAn ) 
for any n G N and each j = 0, 1, 2, and such that 

A n ^A in i^ iv (K 3 ), 

as n — > oo. The existence of such a sequence follows from the density of 
L 1 ^ 3 ,^ 4 ) n i^div( R3 ) in ^div( E3 )' and statement (ii) in Lemma 2.1. For 
each integer n, we know that 

dJp V (A„)(t), a) = / ( PAn * - 3A n ■ a) , 

for any four-potential (t), a) G C^°(1R 3 ,M 4 ). Combining the continuous dif- 
ferentiability of the functional Tpy with statement (i) in Lemma 6.1, we 
obtain, taking the limit n — > oo, 



dJ r PV (A)(D, o) = / (p A t> - J A ■ a) , 

JIR3 V / 



which completes the proof of (ii) in Theorem 2.2. 

Concerning (Hi), recall that the second order differential of F-p\ is equal 

to 

d 2 Jp V (A) = d 2 J- 2 (F) + d 2 J- 4 (A) + d 2 K(A). 

Since Ti is quadratic and J-4 is quartic, estimate (2.27) appears as a conse- 
quence of Formula (2.20), and inequalities (4.19) and (4.57). This completes 
the proof of Theorem 2.2. □ 



References 

[1] V. Bach, J. M. Barbaroux, B. Helffer, and H. Siedentop, On the stability of 
the relativistic electron-positron field, Commun. Math. Phys., 201 (1999), pp. 445-460. 

[2] V. Bach, E. H. Lieb, and J. P. Solovej, Generalized Hartree-Fock theory and the 
Hubbard model, J. Statist. Phys., 76 (1994), pp. 3-89. 

[3] J. D. Bjorken and S. D. Drell, Relativistic quantum fields, McGraw-Hill Book 
Co., New York, 1965. 

[4] P. Chaix and D. Iracane, From quantum electrodynamics to mean field theory: I. 
The Bogoliubov-Dirac-Fock formalism, J. Phys. B, 22 (1989), pp. 3791-3814. 

[5] P. Chaix, D. Iracane, and P.-L. Lions, From quantum electrodynamics to mean 
field theory: II. Variational stability of the vacuum of quantum electrodynamics in 
the mean-field approximation, J. Phys. B, 22 (1989), pp. 3815-3828. 



DIRAC'S VACUUM IN ELECTROMAGNETIC FIELDS 



55 



[6] F. Della Valle, G. Di Domenico, U. Gastaldi, E. Milotti, R. Pengo, G. Ru- 
OSO, AND G. Zavattini, Towards a direct measurement of vacuum magnetic birefrin- 
gence: PVLAS achievements, Optics Communications, 283 (2010), pp. 4194-4198. 
[7] P. A. DlRAC, Theory of electrons and positrons. The Nobel Foundation, 1933. Nobel 
lecture delivered at Stockholm. 

[8] , Discussion of the infinite distribution of electrons in the theory of the positron, 

Proc. Camb. Philos. Soc, 30 (1934), pp. 150-163. 
[9] , Theorie du positron, Solvay report, XXV (1934), pp. 203-212. 

[10] I. Ekeland and R. Temam, Convex analysis and variational problems, vol. 28 of 
Classics in Applied Mathematics, Society for Industrial and Applied Mathematics 
(SIAM), Philadelphia, PA, english ed., 1999. Translated from the French. 

[11] E. Engel, Relativistic Density Functional Theory: Foundations and Basic Formal- 
ism, vol. 'Relativistic Electronic Structure Theory, Part 1. Fundamentals', Elsevier 
(Amsterdam), Schwerdtfeger ed., 2002, ch. 10, pp. 524-624. 

[12] M. J. Esteban, M. Lewin, and E. Sere, Variational methods in relativistic quan- 
tum mechanics, Bull. Amer. Math. Soc. (N.S.), 45 (2008), pp. 535-593. 

[13] H. Euler, tjber die Streuung von Licht an Licht nach der Diracschen Theorie, Ann. 
d. Phys., 418 (1936), pp. 398-448. 

[14] H. Euler and B. Kockel, Uber die Streuung von Licht an Licht nach der Diracschen 
Theorie, Naturwissenschaften, 23 (1935), pp. 246-247. 

[15] F. Finster, A formulation of quantum field theory realizing a sea of interacting Dirac 
particles, Lett. Math. Phys., 97 (2011), pp. 165-183. 

[16] W. M. Frank and A. L. Licht, The nature of perturbation expansion in regularized 
field theories, Nuovo Cimento (10), 31 (1964), pp. 682-687. 

[17] W. Furry, A symmetry theorem in the positron theory, Phys. Rev., 51 (1937), p. 125. 

[18] P. Gravejat, M. Lewin, and E. Sere, Ground state and charge renormalization in 
a nonlinear model of relativistic atoms, Commun. Math. Phys., 286 (2009), pp. 179- 
215. 

[19] , Renormalization and asymptotic expansion of Dirac 's polarized vacuum, Com- 
mun. Math. Phys., 306 (2011), pp. 1-33. 

[20] W. Greiner and J. Reinhardt, Quantum electrodynamics, Springer, 2008. 

[21] W. Greiner and S. Schramm, Resource Letter QEDV-1: The QED vacuum, Amer- 
ican Journal of Physics, 76 (2008), pp. 509-518. 

[22] C. Hainzl, M. Lewin, and E. Sere, Existence of a stable polarized vacuum in the 
Bogoliubov-Dirac-Fock approximation, Commun. Math. Phys., 257 (2005), pp. 515- 
562. 

[23] , Self- consistent solution for the polarized vacuum in a no-photon QED model, 

J. Phys. A, 38 (2005), pp. 4483-4499. 

[24] C. Hainzl, M. Lewin, E. Sere, and J. P. Solovej, A minimization method for 
relativistic electrons in a mean-field approximation of quantum electrodynamics, Phys. 
Rev. A, 76 (2007), p. 052104. 

[25] C. Hainzl, M. Lewin, and J. P. Solovej, The mean-field approximation in quan- 
tum electrodynamics: the no-photon case, Comm. Pure Appl. Math., 60 (2007), 
pp. 546-596. 

[26] C. Hainzl and H. Siedentop, Non-perturbative mass and charge renormalization in 
relativistic no-photon quantum electrodynamics, Commun. Math. Phys., 243 (2003), 
pp. 241-260. 

[27] W. Heisenberg and H. Euler, Folgerungen aus der Diracschen Theorie des 

Positrons, Zeitschrift fur Physik, 98 (1936), pp. 714-732. 
[28] R. Karplus and M. Neuman, Non-linear interactions between electromagnetic 

fields, Phys. Rev., 80 (1950), pp. 380-385. 
[29] T. Kato, Perturbation theory for linear operators, Springer, second ed., 1995. 
[30] G. Leibbrandt, Introduction to the technique of dimensional regularization, Rev. 

Mod. Phys., 47 (1975), pp. 849-876. 
[31] G. Nenciu and G. Scharf, On regular external fields in quantum electrodynamics, 

Helv. Phys. Acta, 51 (1978), pp. 412-424. 



56 



P. GRAVEJAT, C. HAINZL, M. LEWIN, AND E. SERE 



[32] W. Pauli and F. Villars, On the invariant regularization in relativistic quantum 

theory, Rev. Modern Physics, 21 (1949), pp. 434-444. 
[33] M. Peskin and D. Schroeder, Introduction to quantum field theory, Advanced 

Book Program, Addison- Wesley Pub. Co., 1995. 
[34] M. Reed and B. Simon, Methods of Modern Mathematical Physics. I. Functional 

analysis, Academic Press, 1972. 
[35] , Methods of Modern Mathematical Physics. II. Fourier analysis, self- 

adjointness, Academic Press, New York, 1975. 
[36] J. Schwinger, On gauge invariance and vacuum polarization, Phys. Rev. (2), 82 

(1951), pp. 664-679. 

[37] E. Seiler and B. Simon, Bounds in the Yukawa^ quantum field theory: upper bound 
on the pressure, Hamiltonian bound and linear lower bound, Commun. Math. Phys., 
45 (1975), pp. 99-114. 

[38] R. Serber, Linear modifications in the Maxwell field equations, Phys. Rev. (2), 48 
(1935), pp. 49-54. 

[39] B. Simon, Trace ideals and their applications, vol. 35 of London Mathematical Society 
Lecture Note Series, Cambridge University Press, Cambridge, 1979. 

[40] D. A. Slavnov, A generalized Pauli-Villars regularization, Teoret. Mat. Fiz., 17 
(1973), pp. 342-358. 

[41] , Generalized Pauli- Villars regularization in the presence of zero-mass particles, 

Teoret. Mat. Fiz., 19 (1974), pp. 3-13. 
[42] , A direct proof of the correctness of the generalized Pauli- Villars regularization, 

Teoret. Mat. Fiz., 25 (1975), pp. 43-48. 
[43] B. Thaller, The Dirac equation, Texts and Monographs in Physics, Springer- Verlag, 

Berlin, 1992. 

[44] E. Uehling, Polarization effects in the positron theory, Phys. Rev. (2), 48 (1935), 
pp. 55-63. 

[45] J. C. Ward, An identity in quantum electrodynamics, Phys. Rev., 78 (1950), pp. 182- 
182. 

[46] V. Weisskopf, Uber die Elektrodynamik des Vakuums auf Grund der Quantentheorie 
des Elektrons, Math.-Fys. Medd., Danske Vid. Selsk., 16 (1936), pp. 1-39. 

Centre de Mathematiques Laurent Schwartz (UMR 7640), Ecole Polytech- 
nique, F-91128 Palaiseau Cedex, France. 

E-mail address: gravejat@math.polytechnique.fr 

Mathematisches Institut, Auf der Morgenstelle 10, D-72076 Tubingen, Ger- 
many. 

E-mail address: christian.hainzl@uni-tuebingen.de 

CNRS & Laboratoire de Mathematiques (UMR 8088), Universite de Cergy- 
Pontoise, F-95000 Cergy-Pontoise, France. 
E-mail address: mathieu.lewin@math.cnrs.fr 

Ceremade (UMR 7534), Universite Paris-Dauphine, Place du Marechal de 
Lattre de Tassigny, F-75775 Paris Cedex 16, France. 
E-mail address: sere@ceremade.dauphine.fr 



